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EXECUTIVE  SUMMARY 


Alternative  refrigerants  are  being  developed  by  industry  to  prevent  the  further  destruction  of 
stratospheric  ozone  by  chloro fluorocarbons  (CFCs),  which  had  been  the  working  fluids  of  choice  for 
many  air-conditioning  and  refrigeration  machines.  Hydrofluorocarbons  (HFCs)  are  one  class  of 
compounds  that  are  being  pursued  as  replacements  because  their  ozone  depletion  potential  is  zero.  In 
general,  the  exchange  of  fluorine  atoms  on  an  HFC  molecule  with  hydrogen  atoms  decreases  its 
atmospheric  lifetime,  and  it  may  also  increase  the  efficiency  of  the  working  fluid.  Both  of  these  effects 
are  highly  desirable  from  environmental  considerations  since  they  act  to  mitigate  global  warming. 

Unfortunately,  more  hydrogen  on  an  HFC  is  usually  associated  with  an  increase  in 
flammability.  An  accepted  method  for  determining  the  flammability  limits  of  gaseous  fuels  is  ASTM 
Standard  E 681.  The  minimum  and  maximum  concentrations  of  the  fuel  in  air  for  flame  propagation 
are  based  upon  the  observed  ignition  and  growth  of  a flame  in  a vessel  filled  with  a quiescent  fuel/air 
mixture.  A clear  distinction  is  sought  between  a non-propagating  flicker  and  a flame  which  has  enough 
horizontal  propagation  to  be  hazardous.  When  applied  to  hydrocarbons,  these  tests  give  well-defined 
results.  Weak  fuels  like  many  HFCs  have  a great  sensitivity  to  the  test  conditions  and  provide 
ambiguous  limits.  The  ignition  source,  vessel  geometry  and  operator  subjectivity  contribute  to  this 
ambiguity. 

Many  of  the  difficulties  associated  with  the  ASTM  apparatus  are  not  present  in  a planar,  twin- 
flame,  counter-flow  burner.  With  this  burner,  steady  flames  first  are  established  under  favorable 
conditions,  and  then  the  fraction  of  fuel  is  diminished  in  small  increments  until  the  flame  is 
extinguished.  The  tests  are  repeated  at  progressively  lower  flow  rates,  which  are  inherently  more 
capable  of  sustaining  combustion.  By  plotting  the  fraction  of  fuel  at  extinction  versus  the  flow  rate  and 
extrapolating  to  an  experimentally  unattainable  zero-flow  condition,  an  unambiguous  limit  of 
flammability  can  be  attained.  Unlike  the  ASTM  apparatus,  the  counter-flow  burner  method  entirely 
avoids  issues  surrounding  the  design  of  an  ignition  mechanism,  it  minimizes  heat  loss  and  wall  effects, 
and  also  it  is  amenable  to  computational  analysis. 

The  Air-Conditioning  and  Refrigeration  Technology  Institute  began  a program  at  the  Building 
and  Fire  Research  Laboratory  of  National  Institute  of  Standards  and  Technology  in  October,  1994  to 
examine  the  potential  for  the  counter-flow  burner  to  accurately  measure  flammability  limits  of 
refrigerants.  In  phases  I and  II  of  this  research  (reported  earlier),  the  feasibility  of  using  the  counter- 
flow  burner  to  define  the  lower  flammability  limit  (LFL)  of  R-32  (CH2F2)  and  the  minimum  ratio  of  the 
moles  of  R-125  (C2HF5)  to  moles  of  R-32  to  render  an  R-125/R-32  mixture  non-flammable  (i.e.,  the 
critical  flammability  ratio,  CFR)  were  demonstrated.  Theoretical  calculations  were  made  of  the 
structure  of  R-32  flames  in  air  and  in  mixtures  with  R-125  to  assist  in  interpretation  of  the  results.  The 
objectives  for  phase  III  of  the  study  were  (i)  to  determine  the  sensitivity  of  the  measured  LFL  to 
changes  in  burner  design,  (ii)  to  modify  the  burner  to  operate  at  elevated  temperatures  and  with  liquid 
refrigerants  (R-245ca,  or  C3H3F5),  and  (iii)  to  design  a simple  burner  that  could  be  used  by  industry 
and  other  research  laboratories  for  estimating  the  flammability  of  new  refrigerants. 

The  most  significant  design  variable  that  affects  the  LFL  is  the  burner  diameter.  The 
flammability  limits  included  in  this  report  were  obtained  using  a 12  mm  ± 0.2  mm  diameter  contoured 
nozzle  burner,  with  12  mm  ±0.2  mm  spacing.  A 17  % increase  in  burner  diameter  makes  it  difficult 
to  stabilize  the  R-32  flame,  while  a 17  % decrease  in  diameter  results  in  a 20  % decrease  in  the  LFL. 
The  repeatability  of  the  data  for  a fixed  burner  geometry  is  high,  with  results  varying  less  than  ±3  % 
of  the  LFL  for  R-32  and  the  CFR  for  R-125/R-32.  The  variability  of  the  upper  flammability  limit  for 
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R-32  is  about  ± 4 %.  The  greatest  variation  (approximately  ± 10  %)  is  in  the  conditions  which  lead  to 
a flammable  mixture  of  R-245ca  and  air. 

The  following  flammability  limits  (in  terms  of  mole  fraction)  and  critical  flammability  ratios 
(moles  R-125  to  R-32)  were  found  during  the  course  of  this  research,  all  with  air  at  a pressure  of  98 
kPa  ± 2 kPa: 

R-32.  lower  limits  R-32.  upper  limit 

0.14  ± 0.006*  (dry,  30  °C)  0.27  ± 0.01  (8.5  mgwater/gdiyair,  100  °C) 

0.13  ± 0.004  (8.5  mgwater/gdryair,  100  °C) 

R-125 /R-32.  CFR  R-245ca.  flammable  range 

17.6/82.4  to  18.4/81.6  (dry,  30  °C)  0.10  to  0.12  (8.5  mg^/g^  100  °C) 

21.6/78.4  to  22.4/77.6  (8.5  mg^Jg^ air,  100  °C)  nonflammable  (8.5  mg^Jg^^,  50  °C) 

nonflammable  (dry,  100  °C) 

The  absolute  accuracy  of  the  measured  limits  is  three  to  four  times  more  uncertain  than  the 
repeatability,  due  to  the  sensitivity  to  burner  diameter.  Additional  experiments  with  finer  control  on 
the  burner  spacing  and  with  more  flexible  velocity  profiles  are  needed  to  unravel  the  relation  between 
the  LFL  and  the  flow  field/heat  transfer  in  the  burner.  A two-dimensional  numerical  model  of  the 
burner,  with  the  ability  to  vary  the  geometry  and  inlet  conditions,  is  also  necessary  to  properly  interpret 
the  results.  Unlike  for  the  ASTM  apparatus  (for  which  the  absolute  accuracy  has  never  been 
established  for  hydrofluorocarbons),  a reliable  model  can  be  developed  for  the  counter-flow  burner 
with  a reasonable  amount  of  effort. 

The  research  burner  operates  easily  at  100  °C  with  a controlled  humidity.  It  is  estimated  to 
take  less  than  eight  hours  and  500  mg  of  each  refrigerant  to  determine  the  LFL  and  CFR  of  materials 
with  properties  similar  to  R-32  and  R-125.  About  three  times  the  material  and  time  is  expected  for  low 
flammability  liquid  refrigerants  similar  to  R-245ca.  A design  is  proposed  for  an  industrial  burner 
which  is  simpler  to  fabricate  and  maintain. 

This  report  reviews  the  past  work  done  on  premixed,  counter-flowing  flames,  describes  the 
current  counter-flow  burner  facility  and  operating  procedures,  presents  the  experimental  results  with 
the  analysis  that  yields  the  above  flammability  limits,  and  recommends  further  activities  that  could  lead 
to  a science-based  methodology  for  assessing  the  risk  of  fire  from  refrigeration  machine  working  fluids. 


[Note:  This  document  in  its  entirety  has  also  been  published  as  DOE/CE/23810-98,  under  a contract 
from  the  Air-Conditioning  & Refrigeration  Technology  Institute,  ARTI  MCLR  Project  Number  660- 
52401,  Mod.  3.] 


* The  uncertainty  reported  in  DOE/CE/23810-98  was  improperly  listed  as  ± 0.004. 
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LEAN  FLAMMABILITY  LIMIT  AS 
A FUNDAMENTAL  REFRIGERANT  PROPERTY 
Phase  III 


INTRODUCTION 

Background 

Chlorofluorocarbons  (CFCs)  have  been  phased  out  of  production  for  use  as  refrigerants  because  of 
chlorine's  destructive  effect  on  stratospheric  ozone.  Alternative  chemical  compounds  are  being 
adopted  in  which  the  chlorine  atoms  are  replaced  with  fluorine  and/or  hydrogen  atoms;  however, 
replacing  halogen  atoms  in  the  molecule  with  hydrogen  atoms  can  transform  a nonflammable  CFC  into 
a potentially  flammable  material.  An  example  is  the  nonflammable  R-12  (CC12F2)  and  flammable  R-32 
(CH2F2).  The  efficiency  and  flammability  of  R-32  and  other  alternative  refrigerants  must  be  carefully 
balanced  in  mixtures  to  provide  safe  (nonflammable)  and  environmentally  friendly  (energy  efficient  and 
reduced  or  zero  ozone  depleting  potential)  refrigerants. 

The  safety  of  refrigerants  used  by  the  air-conditioning  and  refrigeration  industry  is  classified 
according  to  rules  set  in  ANSI/ASHRAE  Standard  34-1992  [1],  which  specifically  refers  to  ASTM 
Standard  E 681-1994  [2]  for  determining  the  minimum  concentration  of  a gas  in  air  necessary  to  sustain 
a flame.  There  are  three  classes  of  flammability.  A Class  1 refrigerant  shows  no  ability  to  propagate  a 
flame  when  mixed  with  air  at  101  kPa  and  18  °C.  A Class  2 refrigerant  can  propagate  a flame  when 
mixed  with  air  at  101  kPa  and  21  °C,  but  only  at  concentrations  greater  than  0.10  kg/m3.  It  also  must 
have  a heat  of  combustion  less  than  19  MJ/kg.  A Class  3 flammability  rating  applies  to  a refrigerant 
with  a higher  heat  of  combustion,  or  with  the  ability  to  propagate  a flame  when  mixed  in  concentrations 
less  than  0.10  kg/m3.  The  determination  of  the  minimum  concentration  for  flame  propagation,  or  lean 
flammability  limit  (LFL),  is  based  upon  the  observed  ignition  and  growth  of  a flame  in  a 5 L to  12  L 
vessel  filled  with  a quiescent  fuel/air  mixture  at  a specified  uniform  temperature  and  pressure.  The 
standard  test  attempts  to  draw  a clear  distinction  between  a mixture  which  creates  a non-propagating 
flicker  and  a flame  which  has  enough  horizontal  propagation  to  be  hazardous.  When  applied  to  fuels 
like  methane  or  propane,  these  tests  give  well-defined  results.  Weak  fuels  like  R-32,  however,  have  a 
greater  sensitivity  to  the  test  conditions  and  provide  ambiguous  limits.  The  ignition  source  (spark, 
match-head,  heated  wire)  and  energy  level,  the  complex  geometry  of  the  flame,  and  wall  effects  all 
contribute  to  this  ambiguity. 

Many  of  the  difficulties  associated  with  the  ASTM  measurements  of  the  LFL  are  not  present  in 
the  approach  suggested  by  Law,  et  al.[3]  They  used  a planar,  twin-flame,  counter- flow  arrangement  to 
determine  the  volume  flow  rate  at  which  burning  CH4/air  and  C3H8/air  mixtures  are  extinguished.  By 
repeating  the  experiments  at  diminishing  fuel  concentrations,  it  is  possible  to  plot  the  concentration 
versus  the  flow  rate,  and  to  extrapolate  the  results  to  identify  the  concentration  corresponding  to  an 
experimentally  unattainable  zero-flow  condition,  which  is  an  equivalent  definition  of  the  LFL.  Unlike 
the  ASTM  apparatus,  the  counter-flow  burner  method  entirely  avoids  issues  surrounding  the  design  of 
an  ignition  mechanism,  it  minimizes  heat  loss  and  wall  effects,  and  also  it  is  amenable  to  computational 
analysis. 

In  the  previous  and  current  phases  of  the  work  described  in  this  report,  a twin-flame  counter- 
flow burner,  similar  to  that  used  in  reference  [3],  was  selected  to  examine  the  behavior  of 
refrigerant/air  mixtures.  Figure  1 is  a photograph  of  flames  stabilized  about  the  stagnation  region  of 
two  identical,  vertically-aligned,  counter-flowing  jets  of  premixed  R-32  and  air  under  different  sets  of 


Stoichiometric  Flame  Velocity  = 42  cm/s 


Rich  Flame  Velocity  = 24  cm/s 


Lean  Flame  Velocity  =18  cm/s 


Figure  1.  Photographs  of  opposed-flow  burner  with  twin  flames  from  an  R-32/air  mixture, 

showing  impact  of  varying  flow  velocity  and  equivalence  ratio 
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conditions.  The  flames  can  be  extinguished  in  one  of  two  ways:  by  reducing  the  rate  of  chemical 
reaction  to  increase  the  time  required  for  combustion  to  occur,  or  by  increasing  the  flow  velocity  to 
decrease  the  residence  time  in  the  reaction  zone.  When  the  spacing  between  the  burner  jets  and  the  exit 
velocity  is  fixed,  decreasing  the  concentration  of  fuel  reduces  the  rate  of  chemical  reaction.  At  the 
limiting  lean  mixture,  the  reaction  rate  is  slowed  to  the  point  that  the  residence  time  in  the  flame  (which 
is  maintained  about  constant  since  the  velocity  and  spacing  are  fixed)  is  insufficient  for  complete 
combustion  to  occur,  leading  to  flame  extinction.  On  the  other  hand,  when  a fixed  concentration  of  the 
fuel  is  maintained,  increasing  the  jet  velocity  forces  the  flames  towards  the  stagnation  plane  lying 
equidistant  between  the  two  burner  jet  outlets.  At  a sufficiently  high  velocity  the  time  through  the 
reaction  zone  becomes  so  short  that  the  reactants  pass  through  faster  than  they  can  bum,  decreasing  the 
combustion  efficiency  to  a point  that  not  enough  heat  is  released  to  propagate  the  flame.  Again, 
extinction  follows.  The  residence  time  in  the  burner  can  also  be  controlled  by  changing  the  distance 
between  the  jets.  A lesser  distance  for  a fixed  jet  outlet  velocity  produces  a higher  velocity  at  the 
flame  front,  leading  to  a decrease  in  time  available  for  the  combustion  reaction  to  occur. 

The  effect  of  jet  velocity  and  spacing  can  be  combined  into  a single  parameter  called  the 
global  stretch  rate,  Kg,  defined  as  the  average  velocity  at  the  exit  of  the  burner  jet  divided  by  the 
distance  between  the  exit  plane  and  the  stagnation  plane  (i.e.,  half  the  nozzle  separation).  The 
concentration  (or  mole  fraction)  of  refrigerant  at  extinction  can  be  plotted  against  diminishing  values  of 
Kg,  and  linear  extrapolation  used  to  determine  the  minimum  concentration  of  fuel  required  to 
propagate  a flame  under  the  most  favorable  (or  least  stretched,  i.e.,  Kg  = 0 s"1)  flow  conditions.  This 
is  the  LFL  for  the  refrigerant  for  the  given  initial  conditions  (ambient  temperature,  pressure  and 
relative  humidity). 

The  extrapolation  to  a zero-stretch  condition  approach  had  been  used  prior  to  the  current  study 
to  determine  the  flammability  limits  of  highly  flammable  fuels  such  as  hydrocarbons.  The  LFL  of 
slightly  flammable  fuels,  including  many  proposed  refrigerants,  previously  had  not  been  measured  in  a 
flowing  arrangement.  The  feasibility  of  using  an  opposed-flow  burner  to  measure  the  limiting  lean 
mixture  for  a weakly  flammable  refrigerant  gas  in  air  at  standard  conditions  was  evaluated  in  Phase  I of 
this  project  [4],  Phase  I additionally  proved  the  feasibility  of  using  an  opposed-flow  burner 
arrangement  for  determining  the  critical  flammability  ratio  (CFR)  of  R-125/32  mixtures,  where  the 
CFR  is  the  minimum  mole  fraction  necessary  of  the  non-flammable  component  to  render  a binary 
refrigerant  mixture  non-flammable. 

Phase  II  of  the  refrigerant  flammability  project  [5]  was  devoted  to  designing  and  evaluating  an 
improved  version  of  the  concept  developed  in  Phase  I.  A new  opposed-flow  burner  test  facility  was 
built  to  provide  repeatable,  precisely  described  flames  and  extinction  conditions.  Analysis  of  the 
experiments  indicated  that  it  is  possible  to  determine  the  lower  flammability  limit  of  R-32  repeatable 
to  better  than  1 % of  the  limit  value.  The  absolute  uncertainty  reported,  5 %,  includes  a conservative 
assessment  of  the  impact  of  the  uncertainty  in  the  flow  system  calibration. 

Objectives 

The  ultimate  goal  of  this  research  program  is  to  provide  industry  with  an  accurate,  repeatable  method 
for  measuring  flammability  limits.  Phases  I and  II  produced  the  design  of  a measurement  method  that 
is  a radical  departure  from  ASTM  E 681-1994,  utilizing  an  opposed-flow  burner.  Experimental  data 
taken  in  the  earlier  phases  demonstrated  that  the  lower  flammability  limit  of  R-32/air  mixtures  and  the 
critical  flammability  ratio  of  R-32/R- 125/air  mixtures  could  be  determined  accurately,  resulting  in  data 
consistent  with  those  measured  with  the  ASTM  E 681-1994  apparatus,  and  with  a high  degree  of 
repeatability.  The  following  are  significant  questions  which  remained  about  the  opposed-flow  burner 
apparatus,  the  answers  to  which  are  the  objectives  of  Phase  III: 
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• Is  the  absolute  value  of  the  LFL  sensitive  to  small  changes  in  burner  design? 

• Can  the  burner  operate  and  produce  flame  limits  and  critical  flammability  ratios  at  elevated 
temperatures  [e.g..  Elevated  Temperature  Flammability  Limit  at  100  °C  (ETFL100)  and 
CFR  @ 100  °C],  or  with  marginally  flammable  refrigerants  (i.e.,  ASHRAE  class  2 
refrigerants  and  class  1/2  refrigerant  blends)? 

• How  might  one  simplify  the  burner  design  and  economize  its  operation  for  use  by  industry, 
and  what  are  the  tradeoffs  between  designs? 

The  research  has  been  organized  around  three  distinct  tasks  to  answer  these  questions. 

Sensitivity  of  LFL  to  burner  design:  The  theory  states  that  the  extinction  stretch  rate  can  be 
approached  either  by  increasing  the  velocity  for  a fixed  nozzle  spacing,  or  by  decreasing  the  spacing 
for  a fixed  exit  velocity.  The  first  approach  was  used  in  Phase  II,  with  the  nozzle  spacing  fixed  at  12 
mm.  In  Phase  III,  the  burner  spacing  has  been  set  at  different  values  to  determine  the  effect  on  the 
zero-stretch  limit.  Ideally,  the  absolute  value  of  the  LFL  should  not  change,  but  associated 
perturbations  in  buoyancy  and  heat  transfer  can  influence  the  results.  Additional  experiments  have 
been  conducted  with  the  diameter  of  the  burner  changed,  and  the  LFL  measurement  repeated. 

The  shape  of  the  nozzle  determines  the  velocity  profile.  This  geometric  parameter  also  has 
been  examined,  with  the  converging  nozzle  replaced  with  a straight  tube. 

Burner  operation  at  elevated  temperature  and  with  less  flammable  fuel:  Current  industry  safety 
standards  require  the  measurement  of  elevated  temperature  flame  limits  (ETFL)  at  100  0 C.  The 
opposed  flow  burner  has  been  modified  to  operate  at  inlet  gas  temperatures  up  to  this  level.  The  air 
and  refrigerant  mixture  are  heated  in  the  burner,  the  lower  portion  with  heating  tapes  and  the  upper 
portion  with  water  previously  passed  through  a small  boiler  to  elevate  the  temperature  above  100  °C. 
Thermocouples  located  in  the  burner  monitor  the  outlet  temperature.  The  operation  has  been  evaluated 
at  atmospheric  pressure  with  R-32/humid  air  mixtures  at  100  °C. 

Most  measurements  to  date  in  the  NIST  burner  have  had  R-32  as  one  of  the  components  in  the 
fuel/air  mixture.  While  R-32  is  an  order  of  magnitude  less  flammable  than  a hydrocarbon  like 
propane,  it  is  not  the  least  flammable  "slightly-flammable  refrigerant"  being  considered  by  the  air- 
conditioning  and  refrigeration  industry.  Operation  of  the  opposed-flow  burner  presupposes  that,  if  a 
refrigerant  is  reactive  enough,  a flame  can  be  stabilized  under  close-to-stoichiometric  conditions. 

When  this  is  not  the  case,  a different  operating  procedure  is  required.  The  refrigerant  R-245ca  has 
been  reported  to  have  a finite,  but  small,  flammability  range.  However,  a steady  flame  cannot  be 
obtained  when  R-245ca  is  the  sole  fuel.  R-32  has  been  mixed  in  to  increase  the  ease  of  ignition.  By 
decreasing  the  concentration  of  R-32  for  a fixed  flow  of  R-245ca,  extinction  occurs  for  a particular 
stretch-rate/equivalence  ratio  condition.  This  procedure  is  repeated  and  the  results  extrapolated  to 
identify  the  minimum  amount  of  R-32  required  to  have  a stable  flame  under  ideal  (i.e.,  zero  stretch  and 
close-to-stoichiometric  mixture)  conditions.  Measurements  have  been  made  of  the  ETFLioo  of  R-245ca 
with  air  containing  0.0086  g water  vapor  per  g dry  air  (i.e.,  50  % relative  humidity  at  23  °C). 
Measurements  of  the  CFR  of  an  R- 125/32  blend  at  ambient  temperature  and  100  °C  with  air  containing 
0.0086  g water  vapor  per  g dry  air  are  also  included  in  this  report. 

Simplifying  burner  design  and  operation:  Great  care  has  gone  into  the  design  of  the  NIST 
facility.  The  nozzles  are  specifically  contoured  to  produce  a flat  velocity  profile  at  the  exit.  The 


4 


Table  1.  Thermophysical  properties  of  refrigerants  examined 


Name 

Formula 

Molec. 

Wt. 

g/mole 

aBoiling 

Pt.  @ 
101.3  kPa 

aSat.  Press. 

@ 

25  °C 

bEnthalpy  of 
Combustion 
MJ/kgfuel 

aViscosity, 
kg/m-s 
(25  °C) 

aThermal 
cond. , 
W/m-°C 
(25  °C) 

aSpecific 
Heat  (Cp), 
kJ/kg-°C 
(25  °C) 

R-32 

ch2f2 

52.0 

U 

o 

r- 

t-H 

in 

i 

1691  kPa 

-9.35 

12.18 

0.01407 

0.8430 

R-125 

chf2cf3 

120 

-48.1  °C 

1377  kPa 

-3.92c to 
-4.12d 

13.11 

0.01394 

0.7964 

R-245ca 

chf2- 

cf7ch,f 

134 

25.2  °C 

100.5  kPa 

-8.79ce to 

-9  72d-e 

513.6  (liq) 
10.4  (gas) 

0.09012  (1) 
0.0124(g) 

1.304  (1) 
0.9185  (g) 

a Huber  et  al.  [36] 

bBurgess  et  al.  [7] 

ccof2 

in  products 

^o  COF2 

in  products 

Estimate 

burner  is  designed  to  operate  over  a wide  range  of  inlet  velocities  and  stretch  rates.  Computer  operated 
mass  flow  controllers  are  used  to  maintain  and  monitor  the  flow.  While  these  procedures  are  necessary 
for  a NIST  standard  facility,  a number  of  constraints  on  the  design  and  operating  conditions  may  be 
relaxed  without  compromising  too  much  the  accuracy  of  the  results.  For  example,  replacing  the 
contoured  nozzles  with  straight  pipe  sections  with  the  same  internal  diameter  (12  mm)  greatly  simplifies 
the  construction  and  maintenance  of  the  burner.  It  also  facilitates  the  control  of  temperature  for  the 
refrigerant  air  mixture. 

An  advantage  of  the  current  approach  is  that  some  systematic  errors  in  flow  measurement  are 
canceled  when  the  data  are  reduced.  For  example,  an  error  in  estimating  the  total  volume  flow  to  the 
burner  propagates  in  a linear  fashion  to  the  estimated  extinction  stretch  rate.  But  because  the  LFL  is 
found  by  extrapolating  to  the  zero  stretch  condition,  the  intercept  is  unaffected  by  a change  in  slope 
caused  by  a systematic  error  in  determining  flow  rate.  Similarly,  if  the  same  method  were  used  to 
calibrate  the  air  and  refrigerant  flow  meters,  the  calibration  of  one  can  be  checked  against  the  other, 
minimizing  the  impact  of  systematic  errors  on  the  absolute  value  of  the  volume  fraction  of  refrigerant 
in  the  mixture.  Thus,  a simplified  (and  less  expensive)  flow  control  system  may  be  suitable  for  an 
industrial-grade  unit. 

An  operating  procedure  has  been  provided  along  with  estimates  of  the  total  amount  of  chemical 
required  and  the  time  involved  to  measure  the  LFL  of  a pure  refrigerant  and  the  CFR  of  a binary 
mixture  (at  100  °C,  101  kPa,  and  50%  relative  humidity).  Drawings  of  the  burner  are  provided  for 
those  who  wish  to  build  their  own  test  apparatus.  In  addition,  a design  for  a new  burner  and  flow 
control  system  is  proposed  to  reduce  the  cost  of  the  facility  and  increase  ease  of  operation. 

Thermochemical  Properties  of  Refrigerants 

The  flammability  limits  of  gaseous  mixture  are  determined  by  the  heat  released  during  the  combustion 
reaction,  the  heat  absorbed  by  the  materials  intimately  involved,  and  the  minimum  temperature 
required  to  overcome  an  activation  energy.  Thermodynamic  properties  (see  Table  1)  control  the  first 
and  second  quantities,  while  chemical  kinetics  dictates  the  third.  The  kinetics  are  discussed  later  in  the 
analysis  section;  the  equilibrium  reactions  are  discussed  below. 

The  complete,  stoichiometric  combustion  of  R-32  in  dry  air  is  given  by  the  following 
expression: 


CH2F2  + (02  + 3.76  N2)  - C02  + 2 HF  + 3.76  N2 
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{Rl} 


The  equivalence  ratio,  <t>,  for  an  arbitrary  mixture  is  defined  as  the  actual  R-32/air  mole  ratio  divided 
by  the  stoichiometric  ratio,  1/4.76;  the  mole  fraction,  XR_32,  is  defined  as  the  moles  of  R-32  divided  by 
the  total  moles  in  the  mixture,  which,  for  a flowing  system  like  the  counter-flow  burner  is  equivalent  to 
the  volume  flow  of  R-32  divided  by  the  total  volume  flow.  Thus,  XR_32iStoich=  1/5.76  = 0.1736  is  the 
stoichiometric  mole  fraction  of  R-32  in  dry  air. 

A stoichiometric  oxidation  reaction  can  be  written  for  any  refrigerant,  regardless  of  whether  or 
not  it  is  flammable.  For  R-125  and  R-245ca,  the  high  F/H  ratios  in  the  fuel  lead  to  varying  amounts  of 
two  additional  products  of  combustion:  COF2  (carbonyl  fluoride)  and  CF4  (tetrafluoromethane).  (Note 
that  if  moisture  is  present,  COF2  quickly  converts  to  C02  and  HF.)  Unlike  R-32,  the  assumed  product 
composition  is  sensitive  to  the  final  temperature.  Stoichiometric  reactions  under  dry,  adiabatic 
conditions  lead  to  high  temperatures  that  maximize  the  amount  of  COF2;  i.e., 

CHF2CF3  + (02  + 3.76  N2)  - 0.5  C02  + HF  + COF2  + 0.5  CF4  + 3.76  N2,  {R2} 

and 

CHF2CF2CH2F  + 2.5  (02  + 3.76  N2)  - 2 C02  + 3HF  + COF2  + 9.40  N2.  {R3} 

At  temperatures  well  below  the  adiabatic  condition,  dry  carbonyl  fluoride  decomposes  to  equal  amounts 
of  C02  and  CF4.  All  three  compounds  are  of  equal  concentration  at  equilibrium  when  the  temperature 
is  about  1000  °C.  COF2  levels  become  inconsequential  at  room  temperature,  leading  to  the  following 
alternative  complete  combustion  reactions: 

CHF2CF3  + (02  + 3.76  N2)  - C02  + HF  + CF4  + 3.76  N2,  {R2  } 


and 


CHF2CF2CH2F  + 2.5  (02  + 3.76  N2)  - 2.5  C02  + 3HF  + 0.5  CF4  + 9.40  N2.  {R3  } 

In  either  case,  the  stoichiometric  amount  of  air  remains  fixed. 

The  mole  fraction  of  the  less  (or  non-)  flammable  component  in  a binary  refrigerant  mixture  is 
given  the  symbol  x-  The  stoichiometric,  adiabatic  combustion  of  R-125  or  R-245ca  with  R-32  is 
written  as 

XCHF2CF3  + (1-x)CH2F2  + (02+3.76  N2)  - 

(1-.5x)C02+  (2-x)HF  + xCOF2  + .5xCF,+3.76  N2,  {R4} 


or 


xCHF2CF2CH2F  + (1-x)CH2F2+  (1  + 1.5x)(02+3.76  N2)  - 

(1  + X)C02+  (2+x)HF+  xCOF2  + (3.76+5.64x)N2.  {R5} 

As  with  the  pure  refrigerants,  the  COF2  in  the  products  reduces  to  50  % C02  and  50  % CF4  when  the 
final  temperature  drops  below  about  200  °C. 

The  actual  temperature  and  products  of  combustion  in  a flame  are  controlled  by  the  chemical 
kinetics  of  the  reactants  and  strongly  influenced  by  heat  loss  to  the  surrounding  environment. 
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However,  it  is  convenient  to  first  examine  what  temperature  and  species  concentrations  would  result 
from  idealized  adiabatic,  equilibrium  conditions.  These  can  be  determined  for  different  values  of 
equivalence  ratio  using  the  NASA  chemical  equilibrium  computer  code  [6].  Figure  2 shows  the  results 
for  R-32/dry  air  mixtures,  initially  at  25  °C  and  101  kPa.  The  peak  temperature,  1940  °C,  occurs 
when  $ = 1.05.  At  around  the  same  $ the  C02  reaches  a peak  mole  fraction  of  about  0.13.  The  mole 
fraction  of  HF  at  stoichiometric  conditions  is  close  to  0.30,  and  it  continues  to  increase  with  O.  The 
mole  fraction  of  H20  is  over  two  orders-of-magnitude  less  than  the  mole  fraction  of  C02. 

Figure  3 is  a similar  plot  for  R-125/dry  air  mixtures.  The  maximum  temperature  is  400  °C 
less  than  in  the  R-32  flame,  and  the  mole  fraction  of  COF2  is  second  only  to  HF,  with  CF4  becoming 
more  dominant  in  rich  mixtures.  C02  is  present  but  at  much  lower  values  than  in  R-32/air  mixtures, 
and  water  vapor  exists  only  in  parts  per  million. 

The  adiabatic  equilibrium  temperature  is  directly  dependent  upon  the  enthalpy  of  formation  of 
the  reactants.  For  R-245ca,  no  specific  data  on  the  enthalpy  of  formation  could  be  located;  however, 
by  examining  the  compilation  by  Burgess  et  al.  [7],  it  is  possible  to  estimate  the  contribution  of  the 
different  CHxFy  groups  to  the  enthalpy  of  formation  of  the  parent  compound.  A value  of  AHf  = 7580 
kJ/kg  ± 380  kJ/kg  was  used  for  the  calculations  shown  in  Figure  4.  The  peak  adiabatic  temperature 
varies  between  1775  °C  and  1843  °C,  depending  upon  which  extreme  in  AHf  is  used  in  the  calculation. 
As  with  R-32,  HF  is  the  dominant  product,  followed  by  C02,  and  neither  mole  fraction  is  much 
effected  by  the  choice  of  AHf.  No  appreciable  H20  is  formed.  The  higher  temperatures  favor  smaller 
molecules,  which  explains  the  significant  difference  in  mole  fraction  of  F-atom  with  different  enthalpies 
of  combustion. 


RESEARCH  FACILITY  DESIGN 

The  experimental  facility  consists  of  a number  of  subsystems  as  shown  schematically  in  Figure  5. 
Central  to  the  design  is  the  research  counter-flow  burner.  The  inlet  flow  control  system  maintains  the 
proper  mixture  and  amount  of  incoming  materials  to  the  burner,  including  the  air,  gaseous  and  liquid 
refrigerants,  water  vapor,  and  nitrogen.  The  exhaust  treatment  system  scrubs  the  acid  gases  formed 
during  combustion  from  the  effluent  before  it  is  released  into  a chemical  hood.  The  burner  operation 
and  data  acquisition  are  controlled  with  a personal  computer  and  custom  software.  Most  of  the 
subsystems  were  described  in  the  Phase  II  report  [5];  however,  additional  discussion  is  provided  below 
to  document  changes  in  the  system  or  where  it  is  necessary  to  improve  understanding. 

Counter-flow  Burner 

The  counter-flow  burner,  shown  in  Figure  6,  is  approximately  cylindrical,  about  100  mm  in  diameter 
and  450  mm  high.  It  rests  on  a 12  mm  thick  aluminum  base  plate  300  mm  in  diameter.  There  are 
upper  and  lower  sections  to  the  burner,  both  of  which  are  identical  in  design.  The  sections  are 
connected  by  four  rods  to  center  the  jets  and  maintain  their  axes  parallel.  Ignition  is  provided  manually 
with  a retractable  butane  lighter.  When  the  fuel  concentration  is  within  flammable  limits,  a symmetric 
twin  flame  is  formed  on  either  side  of  the  mid-plane,  as  can  be  seen  in  the  photograph  in  Figure  1. 

The  critical  dimensions  are  the  nozzle  separation  and  the  nozzle  diameter,  both  of  which  are  fixed  at 
12.0  mm  ± 0.2  mm  for  the  majority  of  the  experiments.  Air  premixed  with  the  fuel  enters  the  upper 
and  lower  sections  through  9.5  mm  tubes.  The  flow  is  uniformly  distributed  and  straightened  with  an 
80  mm  inner  diameter  by  100  mm  long  tube  containing  a 35  mm  long  piece  of  4 mm  cell  size 
honeycomb  and  fine  mesh  screens  to  break  up  large  eddies.  The  converging  nozzle  is  designed  from 
two  matched  cubic  contours  following  the  criteria  of  Morel  [8],  with  an  area  contraction  ratio  of  44:1. 
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Figure  2.  Calculated  equilibrium  composition  and  temperature  of  R-32/dry  air  mixtures 
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MOLE  FRACTION 


EQUIVALENCE  RATIO,  O 


Figure  3.  Calculated  equilibrium  composition  and  temperature  of  R-  125/dry  air  mixtures 
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EQUIVALENCE  RATIO,  $ 


Figure  4.  Calculated  equilibrium  composition  and  temperature  of  R-245ca/dry  air  mixtures,  showing 
impact  of  AHf 
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MOLE  FRACTION 


The  purpose  of  the  straightener  and  contraction  is  to  produce  a more  uniform  velocity  profile  with  low 
turbulence  intensity  at  the  exit.  Nitrogen  flows  in  a 1.0  mm  thick  concentric  annulus  to  quench  the 
reactants  as  they  escape  from  the  flame,  to  prevent  the  flame  from  stabilizing  on  the  nozzle  rim,  and  to 
reduce  entrainment  of  air.  To  eliminate  unwanted  air  currents  around  the  burner  and  direct  the  exhaust 
gases  upward,  the  entire  burner  is  enclosed  within  a 300  mm  inside  diameter,  6 mm  thick  Plexiglas 
tube.  Drawings  of  the  critical  burner  components  are  included  in  the  Appendix  F. 

The  burner  flanges  directly  exposed  to  the  flames  are  water-cooled  to  maintain  their  integrity 
and  to  minimize  heat  transfer  back  into  the  nozzle.  Water  flows  through  a copper  tube  coiled  around 
the  upper  chamber  of  the  burner  to  prevent  the  exhaust  gases  from  preheating  the  upper  section. 
Thermocouples  are  located  on  the  centerline  just  upstream  of  the  contraction  nozzles  to  monitor  the 
incoming  mixture  temperatures.  Flammability  measurements  at  temperatures  above  the  ambient  are 
conducted  by  heating  the  air/refrigerant  mixture  and  controlling  the  temperatures  of  the  upper  and 
lower  burner  sections.  Water  is  passed  through  a 600  W electrical  heat  exchanger  to  boost  the 
temperature  at  the  inlet  of  the  upper  cooling  coil  to  about  110  °C.  The  temperature  of  the  reactants  in 
the  lower  section  of  the  burner  is  maintained  at  about  100  °C  by  using  heating  tapes  wrapped  around 
the  main  body.  (Note  that  heating  tapes  cannot  be  used  on  the  upper  section  of  the  burner  because  the 
high  HF  levels  in  the  exhaust  quickly  degrade  fiberglass  insulation.)  Fine  control  of  the  reactant 
temperatures  is  accomplished  by  adjusting  the  flow  of  water  through  the  upper  and  lower  nozzle 
outlets.  (A  side  benefit  of  keeping  the  temperature  of  the  flanges  close  to  100  °C  is  that  it  significantly 
reduces  corrosion  of  the  stainless  steel  by  preventing  HF  from  condensing.) 

Flow  Control  Systems 

Reactant  supply:  The  facility  is  designed  to  accurately  control  the  flow  of  the  combustion  air,  gaseous 
refrigerants  (e.g.,  R-32  and  R-125),  liquid  refrigerants  (e.g.,  R-245ca),  gaseous  hydrocarbons  (e.g., 
CH4  and  C3H8),  liquid  hydrocarbons  (e.g.,  C5H12),  and  nitrogen.  The  gases  are  stored  in  individual 
cylinders  at  their  respective  room  temperature  vapor  pressures.  The  air  is  certified  to  have  a mole 
fraction  of  02  equal  to  0.2110  ± 0.0002,  with  water  and  hydrocarbon  levels  below  10'6.  The  remaining 
components  (N2,  Ar,  and  C02)  are  as  taken  from  the  atmosphere.  The  R-32  is  claimed  by  the 
manufacturer  to  contain  mass  fractions  of  water  and  non-volatile  residues  of  less  than  10'5  and  non- 
condensable contaminates  of  less  than  1.5%  by  volume.  The  R-32  is  stored  as  a liquid  at  room 
temperatyre  and  1.6  MPa.  The  R-125  contains  a minimum  mole  fraction  of  0.995  C2HF5  and  is  stored 
at  1 .4  MPa  and  room  temperature.  Air  and  ethane  are  the  largest  contaminants  in  the  CH4,  but  the 
minimum  purity  is  0.9995  mole  fraction.  Table  2 summarizes  the  chemical  compounds  used  in  the 
study,  their  purity,  and  the  suppliers. 

Gas  flow  controllers:  The  flow  of  the  individual  gas  streams  is  controlled  through  an  MKS,  Inc.,1 
Multi  Gas  Controller  type  647B.  This  instrument  powers,  reads  and  controls  the  mass  flow  controllers 
(MFCs)  listed  in  Table  3.  The  resolution  of  the  gas  controller  was  increased  from  0.1  % to  0.01  % of  the 
full  scale  for  each  controller  by  modifications  to  the  software.  All  flow  controllers  are  calibrated  to±  2 % 
of  value  with  the  gas  used  during  testing.  The  calibration  reference  is  a digital  bubble  meter.  This  flow 
meter  was  tested  using  the  NIST  standard  piston  prover  and  has  been  shown  [9]  to  be  accurate  to  within 
± 1 % over  its  entire  range,  0. 1 L/m  to  25  L/min. 


1 Certain  trade  names  and  company  products  are  mentioned  in  the  text  or  identified  in  an  illustration  in  order  to 
specify  adequately  the  experimental  procedure  and  equipment  used.  In  no  case  does  such  identification  imply 
recommendation  or  endorsement  by  the  National  Institute  of  Standards  and  Technology,  nor  does  it  imply  that  the 
products  are  necessarily  the  best  available  for  the  purpose. 
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Figure  5.  Block  diagram  of  experimental  facility 
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Thermocouple 
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Figure  6.  Overall  view  of  counter- flow  burner 
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Table  2.  Chemical  compounds  used  in  study 


Common  Name 

Chemical  Formula 

Minimum  Purity, 
mole  fraction 

Supplier 

air 

N2,  Ar,  C02 

0.2110  ± 0.0002  02 

Matheson 

R-32 

ch2f2 

0.998 

Allied  Signal 

R-125 

c2hf5 

0.995 

Allied  Signal 

R-245ca 

chf2cf2ch2f 

unspecified 

Allied  Signal 

R-245ca 

chf2cf2ch2f 

0.98 

Lancaster  Synthesis 

methane 

ch4 

0.9995 

Matheson 

propane 

c3h8 

0.95 

Air  Products 

isopentane 

c5h12 

0.95 

Fluka 

nitrogen 

n2 

0.99998 

MG  Industries 

Table  3. 

Mass  flow  controllers 

MFC  # 

Calibrated  Gas 

Approximate  Range, 
L/m 

Uncertainty, 

L/m 

MKS  Model  Number 

1 

Air  (top) 

5.0 

± 0.05 

1159B-05000-SN 

2 

R-32  (top) 

2.5 

± 0.012 

1359C-05000-SN 

3 

R-125  (top) 

0.15 

± 0.0008 

1359C-00500-SN 

4 

R-125  (bottom) 

0.30 

± 0.0015 

1359C-01000-SN 

5 

R-32  (bottom) 

1.0 

± 0.005 

1359C-02000-SN 

6 

Air  (bottom) 

10.0 

± 0.05 

1359C-10000-SN 
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Gas  pressures  at  the  inlet  of  the  controllers  are  regulated  at  140  kPa  ± 3 kPa  to  standardize 
initial  flow  conditions  and  validate  the  gas  calibrations.  (The  R-32  supply  bottle  is  connected  to  a lower 
pressure  expansion  tank  to  minimize  pulsations  in  the  flow  of  gas.)  Teflon  tubing,  6.4  mm  in  diameter, 
is  used  to  transport  the  gases.  The  more  and  less  flammable  refrigerants  are  mixed  just  after  leaving  the 
flow  controllers.  The  air  is  routed  through  a humidity  system  prior  to  reaching  the  burner.  At  the 
burner  inlet,  the  air  and  the  refrigerant  mixture  are  combined  and  delivered  to  the  upper  and  lower 
sections  of  the  burner  through  9.5  mm  stainless  steel  tubes.  Figure  7 is  an  overall  schematic  of  the 
flow  control  subsystem  showing  how  the  gas  supply  is  connected  to  the  burner. 

Humidity  control:  The  initial  design  for  adding  water  vapor  to  the  air  stream  was  to  bubble  the 
outlet  of  the  mass  flow  controller  through  two  glass  bubblers  operating  in  series.  Each  bubbler  was 
filled  with  0.50  L of  de-ionized  water  in  which  1.25  kg  of  Mg(N03)2  salt  had  been  dissolved.  (The 
magnesium  nitrate  was  supplied  by  Fisher  Scientific  with  an  assay  value  greater  than  99.9%.) 
Greenspan  [10]  has  shown  that  at  25  °C  and  101  kPa,  this  salt  solution  equilibrates  at  a relative 
humidity  of  52.9  % ± 0.2  %.  The  bubblers  were  placed  in  a room  temperature  water  bath  and  the 
humid  air  at  the  outlet  passed  through  a particle  trap  to  ensure  no  water  drops  remained  suspended  in 
the  air  stream. 

The  method  described  by  Greenspan  is  for  a closed  system  in  which  equilibrium  easily  can  be 
obtained.  However,  in  the  open  system  used  here,  the  continuous  introduction  of  dry  air  carries  with  it 
a small  but  constant  amount  of  C02.  The  purity  of  the  solution  decreases  as  the  H20  and  C02  react  to 
form  H2C03,  dropping  the  equilibrium  relative  humidity.  During  the  first  several  minutes  the  humidity 
[as  measured  within  2 % using  a fast  response  digital  hygrometer  (Fisherbrand  Instant  Hygrometer 
#11-661-7B)]  drops  by  0.25  % per  minute.  This  decline  continued  throughout  testing,  making  a 
continual  replenishment  of  the  salt  solution  necessary. 

To  reduce  the  variability  in  the  relative  humidity  and  eliminate  the  need  to  refresh  the 
solutions,  an  alternative  design  was  chosen  that  is  more  appropriate  for  a flowing  system.  A relative 
humidity  of  50  % at  23  °C  and  standard  pressure  corresponds  to  a humidity  ratio  of  0.0086  gwate/gdry  air- 
This  same  humidity  ratio  is  reached  when  air  is  saturated  at  12  °C.  Saturation  is  achieved  by  bubbling 
air  through  pure  de-ionized  water  that  is  cooled  in  a refrigerated  bath.  A single  2.5  liter  bubbler  with 
0.6  L of  water  is  used  for  each  air  stream.  The  outlet  flow  of  cooled,  humidified  air  passes  through  a 
Teflon  condenser  tube  9.5  mm  in  diameter  and  about  1.5  m long  maintained  at  the  same  temperature  to 
trap  water  droplets  that  might  be  carried  through  the  bubblers.  A ± 0.7  °C  swing  in  temperature 
causes  the  equilibrium  humidity  ratio  to  vary  ± 0.0004  g^Jg^^.  Thus,  the  humidity  ratio  can  be 
increased  or  decreased  by  increasing  or  decreasing  the  temperature  of  the  bath,  even  though  the  air 
stream  may  not  be  fully  equilibrated  with  moisture.  The  temperature  in  the  bubbler  is  measured  with  a 
mercury  thermometer  to  within  0.2  °C,  and  the  humidity  is  checked  prior  to  each  test  with  the 
hygrometer  at  the  inlet  to  the  burner.  A sketch  of  the  humidifier  is  shown  in  Figure  8. 

Liquid  refrigerant  flow  control:  Liquid  fuels  and  refrigerants  with  boiling  points  near  room 
temperature  can  be  either  pre- vaporized  and  metered  as  a gas  or  pre-cooled  and  metered  as  a liquid. 

The  latter  technique  was  chosen  for  the  R-245ca.  Two  100  ml  stainless  steel  syringes  with  silicon 
rubber  O-ring  seals  are  used  to  store  the  refrigerant.  The  syringes  are  encased  in  copper  cooling  coils 
to  prevent  the  refrigerant  from  boiling.  (Pyrex  glass  syringes  with  Teflon  seals  and  a steel  body  were 
found  to  be  unsatisfactory  because  uneven  thermal  expansion  caused  them  to  leak  profusely  at 
temperatures  below  15  °C.)  A Harvard  Apparatus  model  975  syringe  drive  depresses  the  plungers  to 
deliver  between  0.1  and  2.0  ml/min  to  within  ± 0.01  ml/min.  The  temperature  as  measured  with  a 
thermocouple  on  the  outside  of  the  syringe  is  controlled  to  between  4 °C  and  7 °C  by  circulating  water 
from  an  Endocal  RTE-100  chiller  through  the  copper  coils. 
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Figure  7.  Schematic  of  gas  supply  and  flow  controls  for  flammability  measurements 
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Figure  8.  Components  of  air  humidifier 
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The  liquid  refrigerant  flows  from  the  syringe  in  a 3 mm  stainless  steel  tube  directly  to  a 13  mm 
diameter,  150  mm  long  stainless  steel  tube  filled  with  glass  beads  and  wrapped  with  a 200  W heating 
tape.  The  temperature  inside  the  tube  is  heated  to  approximately  80  °C  to  fully  vaporize  the  liquid. 

The  refrigerant  is  then  transported  into  the  burner  through  6.4  mm  Teflon  tubing,  maintained  above  45 
°C  to  prevent  the  refrigerant  from  condensing  in  the  lines  prior  to  injection  into  the  combustion  air 
flow.  A ball  valve  at  the  outlet  of  the  syringe  allows  the  syringe  to  be  filled  from  the  refrigerant  supply 
tank,  as  shown  in  Figure  9. 

R-245ca  is  currently  available  only  in  limited  research  quantities  at  correspondingly  high 
prices.  Because  of  this,  isopentane,  with  a normal  boiling  point  close  to  R-245ca,  was  chosen  as  a 
surrogate  to  test  that  the  syringe  drive  and  boiler  operated  in  a steady  manner.  Even  though  isopentane 
is  much  more  flammable  than  R-245ca,  the  steadiness  of  the  isopentane  flame  was  a direct  indication 
that  the  R-245ca  could  be  delivered  precisely  and  continuously  during  the  performance  of  a 
flammability  measurement. 

The  syringe  drive  was  also  tested  with  decafluoropentane  (C5H2FI0),  a non-flammable  solvent 
made  by  DuPont  that  has  a normal  boiling  point  of  50  °C.  An  R-32/air  flame  was  first  established  at 
the  desired  stretch  rate  and  equivalence  ratio.  The  syringe  drive  was  then  started  and  the  R-32  flow 
was  slowly  decreased  until  the  flame  was  extinguished  at  the  critical  flammability  ratio. 

Exhaust  gas  cleanup:  Hydrofluoric  acid  makes  up  30  % of  the  exhaust  stream  from 
stoichiometric  combustion  of  R-32/air  mixtures.  Much  of  the  HF  is  removed  using  a water  spray 
scrubber.  The  components  of  the  exhaust  gas  clean-up  system  are  shown  in  Figure  10.  A 
polypropylene,  corrosion  resistant  blower  (Dayton  model  5C089)  draws  the  combustion  products 
through  polypropylene  pipe  (diameters  between  90  mm  and  140  mm)  and  routes  the  cleaned  gases  to 
the  exhaust  vent  in  a chemical  hood.  A damper  at  the  entrance  to  the  exhaust  pipe  meters  dilution  air 
to  reduce  the  temperature  of  the  gases  and  minimize  disturbances  to  the  flame  caused  by  excessive 
exhaust  vacuum.  A second  damper  placed  just  ahead  of  the  blower  provides  additional  control  on  the 
exhaust  flow.  A nozzle  (BETE  80°  spray  angle  fog  nozzle)  is  used  to  spray  laboratory  water  at  a rate 
of  approximately  550  ml/min  at  69  kPa  (gauge)  counter  to  the  upward  traveling  exhaust  gas.  The 
acidic  water  is  collected  at  the  bottom  of  the  pipe  and  drained  by  gravity  into  a vented,  four  liter 
polyethylene  jar.  About  90  % removal  efficiency  can  be  achieved  under  normal  operating  conditions. 
To  achieve  closer  to  complete  removal,  a second  scrubber  in  series  with  the  first  is  recommended. 

After  one  or  two  flame  extinction  measurements,  the  water  spray  is  stopped,  the  jar  is  removed  and  the 
effluent  is  neutralized  with  NaOH  to  a pH  of  7 ± 1 before  being  disposed  of  in  the  laboratory  sink. 
Extreme  care  must  be  exercised  to  safeguard  the  operator  from  contact  with  the  hydrofluoric  acid 
solution.  The  addition  of  NaOH  to  the  collection  jar  prior  to  operation  may  be  a safer  procedure. 

Operating  Procedure  and  Data  Acquisition 

Overall  experimental  control  is  provided  by  a 486,  33  MHz  personal  computer  with  12  MB 
RAM.  National  Instruments  Lab  VIEW  3.1.1,  operating  with  Microsoft  Windows  3.1,  is  used  to 
communicate  with  the  mass  flow  controllers  over  an  IEEE  Standard  488  GPIB  talker/listener  interface, 
A special  modification  to  the  communication  EPROM  of  the  647B  multi  gas  controller  increases  the 
resolution  of  the  communication  from  0.1%  to  0.01  % of  the  full  scale  of  flow  of  each  MFC.  This 
combined  with  the  high  accuracy  of  the  calibrated  MFCs  allows  the  operator  control  over  each  gas 
flowing  into  the  burner  to  better  than  ± 2 % of  the  reading.  A 12  bit  National  Instrument  Data 
Acquisition  Card  (AT-MIO-16E-10)  collects  readings  from  various  thermocouples,  a cold  junction 
terminal  block  and  a barometer.  The  thermocouple  cold  junction  terminal  block  (NI,  SC-2070) 
provides  the  ambient  temperature  reading,  ± 0.5  °C.  The  atmospheric  pressure  near  the  flame  is 
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Figure  9.  Components  of  liquid  refrigerant  injector 
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Induced  Draft  Blower 
Dayton  model  5C089 
Adjustable  75  c&n  max 


Figure  10.  Schematic  of  exhaust  gas  scrubber 
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continuously  measured  (±  15  Pa)  with  a Druck  145  DPI  Digital  Barometer.  The  computer  controls  the 
settings  and  monitors  ambient,  burner,  and  flow  conditions  through  a custom  programmed  virtual 
instrument  (VI)  called  TWINFLAM.  Programming  is  in  visual  C and  is  specifically  for  running  these 
tests.  Given  a desired  flame  condition,  TWINFLAM  calculates  the  appropriate  gas  flow  values, 
corrects  for  ambient  conditions,  and  sends  signals  via  the  647B  to  the  mass  flow  controllers.  Refer  to 
Appendix  C for  more  information  about  the  data  acquisition  program. 

For  a given  burner  geometry  and  fuel/refrigerant,  either  the  equivalence  ratio,  $,  or  the  global 
stretch  rate,  Kg,  can  be  chosen  as  the  independent  experimental  variable.  The  other  of  the  two  then 
becomes  the  dependent  parameter.  Kg  is  defined  as  the  nominal  velocity  at  the  exit  plane  of  the  nozzle 
divided  by  the  distance  from  the  exit  to  the  stagnation  plane  (or  half  the  nozzle  separation).  The  exit 
velocity  is  calculated  from  the  total  volumetric  flow  through  the  nozzle  (corrected  for  the  water  vapor 
added  to  humidify  the  air  and  the  ambient  pressure  and  temperature  at  the  burner  outlet)  divided  by  the 
exit  area.  The  mole  fraction  of  less  flammable  refrigerant  is  calculated  as  its  volume  flow  divided  by 
the  total  flow  of  the  less  and  more  flammable  refrigerant  mixture.  When  testing  pure  fuels,  this  value 
is  set  to  zero. 

Appendix  A lists  the  step-by-step  operating  procedure  to  measure  the  LFL  and  CFR.  Before 
beginning  an  experiment,  the  burner  is  pre-heated  and  measurements  of  the  air  humidity  are  taken.  An 
initial  lighting  condition  is  chosen  that  is  robust  enough  for  easy  ignition,  but  has  a flow  velocity 
greater  than  the  flame  speed  to  prevent  flashback  into  the  burner.  The  required  flows  for  a given  Kg 
and  $ are  determined  by  TWINFLAM  and  converted  to  GPIB  signals  which  are  sent  to  the  MFCs. 

The  flame  is  ignited  with  a butane  lighter  (Olympian  GM3X).  Depending  on  the  conditions,  the 
luminescent  region  is  10  mm  to  20  mm  in  diameter  and  the  gap  between  the  flames  is  about  4 mm  or 
less.  After  the  initial  flows  are  set,  the  program  goes  into  a monitoring  loop.  While  in  this  mode, 
current  temperature,  pressure  and  flows  conditions  are  read  and  the  updated  values  of  Kg  and  <f>  are  re- 
calculated over  an  average  of  seven  readings.  During  the  test,  the  burner  temperature  is  maintained  by 
adjusting  the  speed  of  the  water  flow  surrounding  the  nozzle  outlets. 

Once  the  flame  has  stabilized,  the  settings  can  then  be  changed  to  alter  the  flame  composition 
and  bring  it  closer  to  extinction.  This  is  done  by  initiating  a "data  read"  which  ends  the  monitoring 
loop  and  appends  the  current  flame  conditions  to  a data  file.  One  or  more  of  the  input  parameters  (Kg , 
O,  % less  flammable  component)  may  then  be  changed  and  the  new  MFC  signals  sent  to  move  the 
flame  toward  extinction.  The  flammability  limits  for  a pure  gas  are  found  by  either  increasing  Kg  or 
decreasing  O to  find  the  lean  limit  (increase  0 for  rich  limits).  To  determine  the  critical  flammability 
ratio  of  gas  mixtures,  the  percent  of  less  flammable  refrigerant  is  increased  throughout  the  test.  At 
each  step-change,  the  flames  move  closer  together  and  then  re-stabiiize  until  the  merged  flames 
extinguish.  When  the  flame  has  reached  the  extinction  point,  the  data  at  this  condition  is  written  to  the 
end  of  the  test  file.  The  gases  are  then  shut  off  and  preparations  are  made  for  the  next  test.  An 
estimate  of  the  amount  of  material  and  time  required  to  conduct  the  tests  is  given  in  Appendix  B. 


EXPERIMENTAL  RESULTS 

A series  of  experiments  were  conducted  to  compare  the  LFL  of  R-32  and  the  CFR  of  R- 125/32 
measured  in  humidified  air  at  100  °C,  to  the  previous  measurements  of  these  quantities  in  dry  air  at 
room  temperature.  In  the  process,  a complete  flammability  curve  for  R-32  in  humidified  air  at  100  °C 
as  a function  of  the  flame  stretch  rate  has  been  generated,  permitting  an  estimate  of  the  rich 
flammability  limit  as  well.  Results  of  the  counter-flow  burner  experiments  with  R-245ca  are  also 
discussed.  Prior  to  giving  the  flammability  results,  however,  the  investigation  into  the  sensitivity  of  the 
measurements  to  changes  in  burner  geometry  is  presented. 


21 


Impact  of  Geometric  Variables  and  Operating  Procedure 

The  flammability  limits  determined  with  the  ASTM  apparatus  are  known  to  vary  with  the  size  of  the 
vessel  and  the  ignitor  geometry,  as  well  as  with  the  details  of  the  operating  procedure.  The  counter- 
flow burner  does  not  depend  upon  a particular  ignitor  design,  but  the  spacing  between  the  nozzle 
outlets  and  the  size  and  shape  of  the  nozzle  affect  the  results.  The  sensitivity  of  the  results  to  changes 
in  operating  procedure  has  been  evaluated.  The  base-line  geometry  had  the  contoured  nozzles  with 
outlet  diameters  of  12  mm,  spaced  12  mm  apart.  This  was  the  standard  condition  for  all  the  earlier 
measurements,  as  well. 

Effect  of  nozzle  shape:  The  contoured  nozzle  is  designed  to  produce  a close  to  flat  velocity 
profile  at  the  exit  of  an  unopposed  jet  discharging  into  a large,  quiescent  volume.  Measurements  of  the 
velocity  profile  2 mm  downstream  of  the  nozzle  exits  were  made  in  a non-reacting  flow  using  a TSI 
IFA  300  hot  wire  anemometer  and  signal  processor.  When  the  two  jets  are  aligned  in  opposition  (in 
contrast  to  discharging  individually  into  a large  quiescent  volume)  the  centerline  speed  is  found  to  be 
well  below  the  nominal  value  estimated  from  the  total  volume  flow  and  the  cross-sectional  area.  This 
can  be  seen  in  Figure  1 1 . Towards  the  edge  of  the  jet,  the  flow  accelerates  to  a value  greater  than  the 
nominal  value,  indicated  by  the  dotted  line.  This  general  shape  has  been  observed  by  others  studying 
the  dynamics  of  opposed  flow  jets  [11].  The  upper  and  lower  flows  were  not  independently  controlled. 
Based  upon  this  result,  the  flow  control  system  was  modified  and  all  subsequent  experiments  were  run 
with  independent  control  of  the  upper  and  lower  jet  velocities. 

The  profile  in  a fully-developed  laminar  pipe  flow  is  parabolic,  producing  a centerline  speed 
that  is  twice  the  average.  If  a constant  area  nozzle  is  used,  the  gas  speed  at  the  center  of  the  outlet  will 
be  somewhat  less  than  twice  the  nominal,  depending  upon  the  length  of  the  nozzle.  The  effect  on  the 
the  measured  LFL  for  CH4  (in  dry  air  at  room  temperature)  of  varying  the  nozzle  design  from  a 
contoured  to  a constant-area  geometry  was  studied  with  the  diameter  held  constant  at  12  mm.  The 
spacing  was  14  mm  for  the  contoured  nozzle  and  14.6  mm  for  the  straight  nozzle.  The  straight  nozzle 
design  increased  the  LFL  to  0.051,  which  is  beyond  the  95  % confidence  interval  of  0.049  ± 0.001  for 
the  contoured  nozzles  spaced  12  mm  apart  and  greater  than  the  value  of  0.048  reported  by  Ishizuka  and 
Law  [27]. 

An  experiment  was  run  with  screens  inserted  into  the  straight  nozzles  to  see  if  this  might 
improve  the  results.  While  the  screens  appeared  to  allow  the  burner  to  be  operated  at  a slightly  lower 
stretch  rate,  they  led  to  an  increase  in  extrapolated  LFL  for  CH4/air  to  more  than  0.053.  It  is  uncertain 
why  this  occurred,  but  if  turbulence  were  generated  by  the  screens  the  flame  may  have  been 
destabilized. 

Effect  of  nozzle  spacing:  Two  1 mm  thick  washers  were  placed  on  each  of  the  four  support 
posts  separating  the  upper  and  lower  nozzles  so  that  the  distance  between  the  two  increased  to  14  mm. 
The  extrapolated  lower  flammability  limit  of  CH4  in  dry  air  at  room  temperature  increased  to  a mole 
fraction  of  0.0497,  which  is  within  the  upper  limit  of  the  95  % confidence  interval  determined  when  the 
spacing  was  12  mm.  The  scatter  in  data  is  slightly  higher  at  the  larger  spacing. 

Effect  of  nozzle  diameter:  The  effect  of  nozzle  opening  was  examined  with  three  different 
diameter  (10  mm,  12  mm,  and  14  mm)  straight  nozzles.  A constant  area  design  was  used  (rather  than 
the  contoured  design)  for  ease  of  fabrication.  With  CH4/air  mixtures,  increasing  the  diameter  from  12 
mm  to  14  mm  (with  the  spacing  kept  constant  at  14.6  mm)  has  no  statistically  significant  effect  on  the 
LFL.  For  the  R-32/air  mixture,  no  stable  flame  could  be  maintained  when  the  diameter  was  increased 
from  12  mm  to  14  mm  (at  a spacing  of  14.6  mm),  and  this  did  not  change  when  screens  were  inserted 
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Figure  11.  Outlet  velocity  measured  in  non-reacting,  counter-flow  jets  from  contoured  nozzles  without 
independent  control.  The  maximum  fluctuation  is  indicated  with  the  error  bar.  Dotted  line  represents 
profile  assuming  uniform  average  velocity  from  total  volume  flow. 
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into  the  straight  nozzle  outlets. 

Decreasing  the  nozzle  diameter  to  10  mm  (with  a nozzle  spacing  of  12.7  mm)  led  to 
significantly  more  stable  flames.  The  extrapolated  LFL  for  CH4/air  decreased  to  about  0.046,  and  the 
LFL  for  R-32/air  (dry  and  at  room  temperature)  dropped  dramatically  from  0.14  to  near  0.11.  The 
general  trend  was  unaffected  by  the  presence  or  absence  of  screens  at  the  nozzle  outlet. 

Conclusions  from  parametric  study  of  geometry:  The  important  conclusions  that  are  drawn 
from  the  parametric  study  are  the  following: 

• The  contoured  nozzles  produce  a slightly  lower  value  of  LFL  than  the  constant  area  nozzles. 

• There  is  no  advantage  to  adding  straightening  screens  at  the  outlet  of  the  nozzle. 

• A spacing  of  about  12  mm  gives  consistent  results  for  LFL  for  both  a hydrocarbon  and  a 
hydrofluorocarbon . 

• R-32/air  flames  are  much  more  sensitive  to  burner  diameter  than  CH4/air  flames,  with  diameters  as 
large  as  14  mm  having  a destabilizing  affect,  and  diameters  as  small  as  10  mm  having  a stabilizing 
effect. 

The  decrease  in  measured  LFL  with  nozzle  diameter  is  problematic  since  it  is  this  sort  of  ambiguity  that 
the  counter-flow  burner  is  meant  to  overcome.  Additional  experiments  with  finer  control  on  the  burner 
spacing  and  with  more  flexible  velocity  profiles  are  needed  to  unravel  the  relation  between  the  LFL  and 
the  flow  field/heat  transfer  in  the  burner.  A two-dimensional  numerical  model  of  the  burner,  with  the 
ability  to  vary  the  geometry  and  inlet  conditions,  is  also  necessary  to  properly  interpret  the  results. 
Unlike  for  the  ASTM  apparatus,  a reliable  model  can  be  developed  for  the  counter-flow  burner  with  a 
reasonable  amount  of  effort.  Pending  further  studies,  one  must  attribute  a measured  LFL  to  a specific 
burner  geometry.  The  12  mm  diameter  contoured  nozzle,  with  12  mm  spacing,  is  recommended  for 
yielding  the  most  consistent  and  reliable  values  for  flammability,  with  an  expected  accuracy  for  that 
particular  geometry  of  ± 3 % of  the  measured  mole  fraction.  In  the  absolute  sense,  however,  the 
uncertainty  will  remain  three  to  four  times  higher  until  additional  research  is  conducted. 

Sensitivity  to  operating  procedure:  The  impact  on  the  LFL  of  the  initial  conditions,  the 
approach  to  extinction  and  the  operation  of  the  exhaust  scrubber  were  examined.  Ignition  can  be 
difficult  if  the  initial  mixture  ratio  and  stretch  rate  are  outside  of  a limited  range.  Flash  back  is  only  a 
problem  for  highly  flammable  fuels  like  methane.  It  was  not  found  to  be  a concern  for  R-32.  The  act 
of  removing  the  ignitor  from  the  lit  flame  can  extinguish  the  flame  before  the  measurement  even 
begins.  Likewise,  if  the  draw  on  the  exhaust  is  too  great,  the  flame  can  be  extinguished  when  the  fan  is 
turned  on.  Once  the  flame  has  reach  a steady  condition  with  the  ignitor  removed  and  the  exhaust 
operating  properly,  there  is  no  memory  of  the  ignition  process,  so  that  the  measured  extinction  value  is 
totally  divorced  from  how  ignition  was  accomplished. 

Buoyancy  dominates  the  flow  at  low  inlet  velocities.  A global  stretch  rate  of  30  s'1  appears  to 
be  a practical  lower  limit  for  producing  an  undistorted  flame,  although  some  flames  can  be  stabilized  at 
a stretch  rate  as  low  as  25  s'1 . When  determining  the  LFL  by  extrapolation,  there  is  no  need  to  run 
tests  at  stretch  rates  greater  than  70  s'1. 

The  mole  fraction  of  refrigerant  at  which  a flame  is  extinguished  for  a fixed  global  stretch  rate 
is  increased  by  external  disturbances  to  the  flame.  Excessive  exhaust  flow,  unintended  leaks  into  the 
Plexiglas  burner  chamber,  and  the  flowing  nitrogen  shield  can  stretch  the  flame  beyond  the  nominal 
value.  All  of  these  effects  become  more  important  at  the  lower  stretch  rates.  A systematic  study  was 
performed  with  variable  flow  in  the  nitrogen  shield,  and  it  was  found  that  if  the  average  velocity  of  the 
nitrogen  gas  was  greater  than  30  % of  the  outlet  velocity  of  the  fuel/air  mixture,  then  the  LFL  was 
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measurably  increased.  Experiments  run  with  zero  nitrogen  flow  gave  consistent  average  results,  but 
the  standard  deviation  was  slightly  higher  for  the  lowest  stretch  rates. 

Flammability  Measurements  of  R-32 

Experiments  were  conducted  in  R-32/air  mixtures  maintained  at  100  °C  ± 5 °C.  The  measurements 
were  all  taken  using  the  12  mm  diameter  contoured  nozzles,  spaced  12  mm  apart,  with  the  nitrogen 
shield  gas  flowing  at  a velocity  about  30  % of  the  reactant  jet  velocity.  Experiments  were  run  with  dry 
air  and  with  the  dew  point  maintained  between  11.6  °C  and  13.4  °C.  A dew  point  of  12.0  °C 
corresponds  to  a relative  humidity  of  50  % at  standard  pressure  and  23  °C. 

Figure  12  is  a plot  of  the  extinction  R-32  mole  fraction  as  a function  of  the  global  stretch  rate. 
The  data  points  are  represented  by  the  open  and  solid  circles,  corresponding  to  dry  and  humidified  air, 
respectively.  The  solid  line  is  a linear  fit  through  the  moist  air  data  for  global  stretch  rates  between  30 
s'1  and  70  s'1.  An  extrapolation  of  the  line  to  a zero  stretch  condition  yields  a lean  flammability  limit 
mole  fraction  (ETFL100)  of  0.131  ± .001.  Also  plotted  in  Figure  12  are  the  data  and  linear  fits  of  the 
room  temperature  measurements  taken  from  the  Phase  II  report  [5]  for  dry  air  (open  diamonds,  dashed 
line)  and  air  with  a relative  humidity  at  standard  conditions  of  43  % ± 2 % (solid  diamonds,  dotted 
line).  Both  the  previous  and  current  work  indicate  that  the  variation  in  LFL  (at  a fixed  temperature)  for 
differing  humidity  is  less  than  the  experimental  uncertainty.  However,  the  temperature  can  be  seen  to 
have  a significant  impact  on  the  LFL,  decreasing  the  lean  limit  from  about  0. 14  for  temperatures  under 
35  °C  to  about  0.13  for  temperatures  near  100  °C.  While  the  intercepts  change  with  the  temperature, 
the  curves  remain  almost  parallel.  The  slope  of  the  extinction  mole  fraction  versus  stretch  rate  curves 
lie  between  0.00046  s and  0.00048  s. 

The  equivalence  ratio  that  leads  to  the  maximum  extinction  stretch  rate  changes  with  the 
properties  of  the  fuel.  For  example,  Law  et  al.  [3]  found  that  for  methane/air  mixtures,  the  highest 
extinction  stretch  rate  occurs  for  o = 0.95,  while  for  propane/air  mixtures,  the  most  robust  mixture  is 
associated  with  an  equivalence  ratio  near  1.20.  To  determine  the  mixture  of  R-32  in  air  that  leads  to 
the  most  difficult  flame  to  extinguish,  additional  experiments  were  performed  for  stoichiometric  and 
rich  conditions. 

Figure  13  is  a plot  of  the  extinction  mole  fraction  for  rich  and  lean  R-32  flames  versus  global 
stretch  rate,  with  the  reactants  initially  at  a nominal  temperature  of  100  °C  and  dew  point  of  12°C. 

The  lower  portion  of  the  curve,  for  Kg  less  than  70  s'1,  includes  the  same  data  as  shown  in  Figure  12. 
The  horizontal  line  corresponds  to  o = 1.0  (XR.32=  0.1736).  It  is  evident  that  the  stoichiometric 
condition  is  nowhere  near  the  most  difficult  to  extinguish.  Operating  at  an  R-32  mole  fraction  of  0.202 
increases  the  extinction  stretch  rate  to  156  s'1,  compared  to  102  s'1  for  a stoichiometric  flame.  The  peak 
mole  fraction  corresponds  to  an  equivalence  ratio  of  about  1.20,  which  is  similar  to  the  behavior  of 
propane/air  flames. 

The  structure  of  rich  flames  is  more  complex  than  flames  in  lean  mixtures  because 
recombination  reactions  lead  to  multi-carbon  species  and  soot,  and  preferential  diffusion  of  the  H-atom 
is  enhanced  among  the  field  of  larger  hydrofluorocarbon  molecules.  Even  so,  inspection  of  the  upper 
branch  of  the  flammability  curve  in  Figure  13  suggests  that  a linear  extrapolation  to  a zero-stretch 
condition  can  yield  an  easily  identifiable  upper  limit.  A straight-line  fit  through  the  data  with  XR_32 
>0.24  has  a y-intercept  of  0.271;  excluding  data  with  stretch  rates  greater  than  60  s'1  produces  a 
slightly  lower  value,  0.263.  These  values  are  less  than  most  of  those  reported  in  the  literature,  which 
range  for  ambient  initial  conditions  from  0.269  using  a flame  tube  and  fuse  wire  [12]  to  0.334 
measured  in  the  ASTM  E-ASTM  E 681-1994  apparatus  with  a match  ignitor  [13]. 
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R-32  Mole  Fraction  in  Mixture 


Figure  12.  Extinction  mole  fraction  of  R-3 2/air  mixtures  as  a function  of  flame  stretch,  showing 
impact  of  initial  temperature  and  humidity.  Solid  symbols:  Tdewpt=  12  °C;  open  symbols:  dry  air. 


26 


R-32  Mole  Fraction 


Figure  13.  Full  flammability  curve  for  R-32/air  mixtures  at  100  °C  and  a 50  % relative  humidity 
(at  23  °C). 
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CFR  of  R-125/32 


The  complete  combustion  at  stoichiometric  conditions  of  a mole  fraction  of  R- 125,  Xr-i25>  in  an  R- 
125/32  refrigerant  mixture  is  shown  in  reaction  {R4}  in  the  background  section.  The  total  equivalence 
ratio  (i.e.,  assuming  both  refrigerants  contribute  to  the  fuel)  is  defined  as 

® = 4.76  (V'R.125  + VW/V'*  (1) 

where  V'  is  the  volume  flow  of  the  respective  gases. 

The  critical  flammability  ratio,  CFR,  is  the  value  of  x required  to  render  the  refrigerant/air 
mixture  non-flammable  under  the  most  conservative  (i.e.,  zero  stretch)  conditions.  The  CFR  was 
measured  during  Phase  I of  this  research  [4]  using  an  earlier  version  of  the  counter-flow  burner.  The 
value  of  the  CFR  for  R-125/32  in  dry  air  (at  a temperature  which  was  not  recorded)  was  estimated  to 
be  0.185  ± 0.008.  These  experiments  were  repeated  using  the  redesigned  research  burner. 

Figure  14  is  a plot  of  the  global  stretch  rate  at  flame  extinction  as  a function  of  total 
equivalence  ratio  for  different  mole  fractions  of  R-125  in  the  refrigerant  mixture  at  about  30  °C.  When 
no  R-125  is  present,  the  maximum  occurs  around  o=  1.25.  The  extinction  stretch  rate  is  not  a strong 
function  of  equivalence  ratio  in  this  region  (for  fixed  x),  but  the  peak  in  the  curve  shifts  downward 
from  when  x=0  to  Opeak<  1.15  when  x=  0.120.  The  mole  fraction  of  R-125  is  plotted  in  Figure  15 
(filled  circles)  as  a function  of  the  corresponding  maximum  extinction  stretch  rate  from  Figure  14.  If  a 
straight  line  is  fit  to  all  the  data  (the  dotted  line),  it  intersects  the  ordinate  when  x(0)=  Xo  = 0.149.  If 
the  zero  stretch  condition  is  identified  by  extrapolating  data  restricted  to  less  than  48  s'1  (the  solid  line), 
then  Xo  is  about  0. 170,  which  can  be  taken  as  a conservative  estimate  of  the  CFR  for  R-125/32 
mixtures  in  dry  air  at  30  °C  ± 5°C  and  99.4  kPa  ± 0.8  kPa. 

The  data  from  the  room  temperature  measurements  of  the  CFR  of  R-125/32  showed  little 
sensitivity  to  the  exact  value  of  the  equivalence  ratio,  but  that  the  most  robust  flames  occurred  between 
an  overall  equivalence  ratio  of  1.15  and  1.25.  To  streamline  the  process  of  obtaining  the  CFR  for 
elevated  temperatures  and  humidity,  the  operating  procedure  was  modified  to  permit  Xr-i25  t0  be 
increased  in  small  increments  while  maintaining  the  overall  equivalence  ratio  and  the  global  stretch  rate 
constant.  Three  values  were  chosen  for  $:  1.15,  1.18,  and  1.20.  The  12  mm  diameter  contoured 
nozzles  with  12  mm  spacing  were  used  in  the  burner. 

The  mole  fraction  of  R-125  at  which  the  flame  extinguished  is  shown  in  Figure  16,  plotted 
versus  the  global  stretch  rate.  The  filled  squares  are  for  $ = 1.15,  the  open  circles  for  $ = 1.18,  and 
the  open  triangles  represent  the  data  from  the  O = 1.20  tests.  The  reactant  temperature  for  this  series 
was  set  to  100  °C  ± 5 °C,  and  the  dew  point  was  kept  between  12.2  °C  and  13.2  °C.  Comparing 
Figure  16  to  Figure  15,  the  stabilizing  effect  of  the  higher  temperature  and  humid  conditions  is  easily 
seen.  The  importance  of  temperature  on  flame  stabilization  is  known  from  theoretical  considerations, 
and  was  also  demonstrated  in  the  R-32  LFL  tests.  The  R-32  experiments  showed  no  sensitivity  to 
relative  humidity,  but  it  is  possible  that  the  water  vapor  becomes  more  significant  in  CFR 
measurements  at  the  higher  values  of  Xr-i25>  where  the  H/F  ratio  becomes  small.  A separate  test  series 
is  required  to  ascertain  the  relative  importance  of  humidity  compared  to  the  temperature  for  enhancing 
the  combustion  as  the  mixture  approaches  the  CFR. 

For  stretch  rates  less  than  30  s'1,  the  value  of  the  R-125  extinction  mole  fraction  falls  off 
quickly,  a result  that  was  not  noted  with  the  dry  air,  room  temperature  measurements.  The  flame  was 
observed  to  operate  in  an  erratic  way  at  these  low  values  of  flame  stretch,  probably  due  to  the  high 
buoyancy  forces  and  interactions  with  the  hotter  nozzle  surface.  Additional  disturbances  might  also  be 
associated  with  the  bubbling  action  of  the  humidifier,  which  was  not  present  in  the  room  temperature 
CFR  study. 
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Global  Stretch  Rate, 


Figure  14.  Extinction  stretch  rates  for  fixed  values  of  x as  a function  of  total  equivalence  ratio 
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0.20 


Figure  15.  Values  of  x required  to  render  R-32/125/air  mixtures  non-flammable  as  a function  of  flame 
stretch 
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Mole  Fraction  R-125 


Global  Stretch  Rate,  s'1 


Figure  16.  Mole  fraction  of  R-125  at  extinction  as  function  of  global  stretch  rate  for  reactants  at 
100  ° C and  dew  point  of  12  °C.  The  dashed  curves  are  best  fits  through  data, 
excluding  high  and  low  values  of  Kg,  corresponding  to  three  values  of  O.  The 
solid  line  is  best  fit  without  regard  to  O. 
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Straight-line  extrapolations  to  the  zero  stretch  condition  were  performed  independently  for  each 
of  the  three  values  of  overall  equivalence  ratio.  The  extinction  values  of  XR-125  for  Kg  less  than  about  30 
s'1  and  greater  than  around  60  s'1  were  not  considered  in  arriving  at  the  dotted  lines  shown  in  Figure 
16.  The  values  of  xo  resulting  from  the  o = 1.15  and  $ = 1.20  tests  are  close  to  each  other;  at  an 
equivalence  ratio  of  1.18,  Xo  reaches  its  maximum  value,  0.230.  The  solid  line  in  Figure  16  is  the  best 
fit  when  all  values  of  4>  are  lumped  together,  yielding  an  intercept  of  Xo  = 0.217.  The  most 
conservative  value  for  the  CFR  at  100  °C  and  a 50  % relative  humidity  (corrected  to  23  °C)  is  0.23. 

Flammability  Measurements  of  R-245ca/R-32  Mixtures 

R-245ca  has  a boiling  point  close  to  room  temperature,  which  necessitates  using  an  alternative  fluid 
metering  system  and  operating  procedure.  Liquid  refrigerant  is  drawn  into  two,  100  ml,  cooled, 
stainless  steel  syringes.  The  syringes  are  placed  onto  a motorized  drive  unit  that  is  equipped  with 
cooling  coils  to  maintain  the  refrigerant  in  a liquid  state,  at  a temperature  below  7.5  °C.  The  speed  of 
the  syringe  drive  is  adjustable  in  discrete  increments  only.  Because  R-245ca  is  much  less  flammable 
than  R-32,  it  is  necessary  to  add  a more  flammable  component  to  the  reactant  stream  to  enable  ignition. 

The  procedure  adopted  to  determine  the  flammability  limits  of  R-245ca/air  mixtures  is  to  first 
ignite  the  flame  with  the  syringe  drive  engaged  but  with  a preponderance  of  R-32  present,  at  an  overall 
equivalence  ratio  around  1.3  and  a stretch  rate  between  50  s'1  and  100  s'1.  Once  the  flame  is  stabilized, 
the  mole  fraction  of  R-32  is  reduced,  while  either  the  overall  equivalence  ratio  or  the  stretch  rate  is 
maintained  constant.  The  fraction  of  R-32  in  the  refrigerant  mixture  and  the  stretch  rate  at  the  point  of 
extinction  are  noted,  and  the  test  is  repeated  for  a different  initial  stretch  rate.  By  plotting  the  mole 
fraction  of  R-245ca  in  the  refrigerant  mixture  as  a function  of  the  extinction  stretch  rate,  a graph 
similar  to  Figure  15  is  obtained. 

The  solid  circles  plotted  in  Figure  17  show  the  results  for  an  initial  reactant  temperature  of  100 
°C  and  a room  temperature  relative  humidity  of  50  %.  The  data  exhibit  a sharp  break  at  a mole 
fraction  of  R-245ca  equal  to  about  0.24.  For  smaller  values  of  Xr-245>  the  refrigerant  acts  in  a manner 
similar  to  R-125,  inhibiting  the  R-32/air  reaction  due  to  its  high  molecular  weight  and  lower  reactivity. 
A straight  line  extrapolation  from  these  low  values  of  xr-245  (solid  line  in  Figure  17)  would  suggest  that 
R-245ca  is  non-flammable,  with  a CFR  in  R-32  of  about  0.34.  On  the  other  hand,  the  zero  stretch 
limit  based  upon  xr-245  > 0-24  (or  Kg  < 60  s'1)  leads  to  the  conclusion  that  R-245ca/air  mixtures  can 
sustain  a flame  under  idealized  conditions  since  a straight  line  intersects  the  y-axis  at  a value  of  Xr-245 
greater  than  unity.  The  data  used  for  extrapolation  represent  overall  equivalence  ratios  between  1.26 
and  1.35. 

The  effects  of  initial  reactant  temperature  and  relative  humidity  on  the  extinction  mole  fraction 
of  R-245ca  can  also  be  seen  in  Figure  17.  The  open  diamonds  represent  tests  done  with  the  initial 
temperature  reduced  to  50  °C  and  the  relative  humidity  held  constant  at  50  %,  and  the  squares  with 
dots  show  what  occurs  if  the  air  is  not  humidified  and  the  temperature  is  maintained  at  100  °C.  Both 
reducing  the  temperature  and  the  relative  humidity  lead  to  lower  values  of  extinction  mole  fraction.  The 
limited  data  taken  at  these  conditions  suggest  that  either  change  is  sufficient  to  cause  the  extrapolated 
straight  line  to  intersect  the  y-axis  below  unity  (xr-245,o  < 1.0),  implying  that  R-245ca  is  non-flammable 
under  these  conditions.  If  a single  correlation  is  derived  using  all  data  for  stretch  rates  less  than  60  s'1, 
the  dotted  line  in  Figure  17  is  the  result. 

Figure  18  shows  how  the  extinction  stretch  rate  is  affected  by  the  overall  equivalence  ratio. 

The  data  are  grouped  according  to  the  value  of  Xr-245-  The  scatter  within  each  grouping  is  associated 
with  the  variations  in  the  initial  temperature,  the  relative  humidity  and  the  exact  value  of  Xr-245-  The 
open  hexagon  symbols  represent  the  data  for  pure  R-32  (xr-245  = 0).  The  strong  influence  of  Xr-245  on 
the  extinction  stretch  rate  is  the  most  obvious  feature  of  the  figure.  A sensitivity  of  Kg  to  overall 
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Figure  17.  Mole  fraction  of  R-245ca  mixed  with  R-32  as  a function  of  extinction  stretch  rate. 

Solid  circles:  100  °C,  50  % relative  humidity;  diamonds:  50  °C,  50  % relative 
humidity;  squares:  100  °C,  0%  relative  humidity. 


33 


Figure  18.  Dependence  of  extinction  stretch  rate  on  total  equivalence  ratio  for  ranges  of  Xr245 
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equivalence  ratio  is  not  observed,  although  a slight  shift  in  the  peak  to  higher  Otot  with  increasing  Xr-245 
can  be  imagined. 

Table  4 summarizes  the  intercept  values  (xr-245.o)  and  correlation  coefficients  found  by  trying  to 
separate  out  the  effects  of  humidity,  temperature  and  equivalence  ratio.  Correlation  coefficients  less 
than  0.90  result  when  the  data  are  scattered  or  not  well  correlated  by  a straight  line.  Intercepts  greater 
than  unity  imply  pure  R-245ca/air  mixtures  are  flammable;  extinction  stretch  rates  less  than  zero  imply 
that  the  R-245ca/air  mixtures  are  non-flammable  at  the  conditions  stated.  Additional  measurements  are 
required  to  better  quantify  the  precise  boundaries  of  flammability  for  R-245ca,  but  clearly  they  are 
sensitive  to  the  initial  conditions,  with  equivalence  ratios  around  1.32,  higher  temperatures,  and  moister 
air  increasing  the  chance  of  flammability.  (The  same  conclusions  regarding  equivalence  ratio  and 
humidity  were  reached  by  Smith,  et  al.  [37],  using  the  ASTM  E 681-1994  apparatus.)  It  is  also 
important  to  keep  the  strict  definition  of  "flammable"  in  perspective,  since  under  no  conditions 
obtainable  in  the  laboratory  counter-flow  burner  could  a pure  R-245ca/air  flame  be  stabilized. 


Table  4.  Intercept  values  for  R-245ca/air  flames  found  using  the  tests  conditions  indicated 


Stretch  rate  s1 

Temperature 

°C 

Relative 
Humidity 
@23  °C 

41  total 

Xr-245, 0 

Correlation 

Coefficient 

Kg.o,  s'1 

30-60 

50-  100 

0 % - 50  % 

1.1  - 1.6 

0.92 

0.77 

< 0 

30-60 

100  ± 5 

50  % ± 4 % 

1.26-  1.35 

> 1.0 

0.96 

6.0 

30-60 

50  ± 5 

50  % ± 4 % 

1.27-  1.41 

0.98 

0.91 

< 0 

30-60 

100  ± 5 

0 % 

1.20-  1.36 

0.58 

0.17 

< 0 

30-60 

50-  100 

0 % - 50  % 

1.26  ± 0.02 

0.91 

0.66 

< 0 

30-60 

50-100 

0 % - 50  % 

1.32  ± 0.02 

> 1.0 

0.70 

12 

30-60 

50-  100 

0 % - 50  % 

1.37  ± 0.02 

> 1.0 

0.89 

3.8 
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ANALYSIS  AND  DISCUSSION 


Uncertainty  Analysis 

The  LFL  and  CFR  reported  in  the  previous  sections  are  subject  to  uncertainties  from  several  sources, 
including  errors  in  flow  measurements;  variations  in  temperature,  pressure  and  composition;  the  rate  at 
which  the  extinction  point  is  approached;  changes  in  burner  geometry;  and  non-linear  effects  near  the 
zero  stretch  rate  condition.  Flow  calibration  and  measurement  uncertainties  were  examined  and 
described  in  detail  for  the  R-32/air  LFL  measurements  in  the  Phase  II  report  [5].  An  uncertainty  in 
mole  fraction  of  ± 0.008,  with  a 95  % confidence  interval,  was  estimated  based  upon  the  assumed 
linear  relation  between  the  LFL  and  stretch  rate,  and  the  dependence  of  each  on  the  uncertainty  in 
measured  flows.  In  the  Phase  III  work,  additional  mass  flow  controllers  (MFCs)  have  been  added 
which  are  sized  to  operate  close  to  the  middle  of  their  dynamic  range  during  most  of  the  experiments. 
The  additional  MFCs  also  allow  independent  control  of  the  upper  and  lower  burner  sections,  to  ensure 
that  the  flows  are  balanced  equally. 

The  impact  on  the  uncertainty  of  the  results  caused  by  variations  in  humidity,  unknown 
concentrations  of  trace  species  in  the  reactants,  and  variations  in  inlet  temperature  can  not  be  expressed 
in  a simple  mathematical  expression  because  of  the  complex  relationships  between  these  parameters  and 
the  flame  chemistry.  The  uncertainty  in  dew  point  is  about  ± 1 °C,  which,  based  upon  the 
measurements  with  and  without  any  moisture  added,  is  estimated  to  be  an  insignificant  source  of 
uncertainty  for  the  ETFL100  of  R-32;  the  error  in  the  CFR@100  °C  of  R- 125/32  and  the  ETFL100  of  R- 
245ca  could  be  larger,  but  is  still  likely  to  be  smaller  than  the  uncertainty  from  the  flow.  The  average 
temperature  of  the  reactants  as  they  enter  the  burner  nozzle  varies  less  than  10  °C,  and  the  maximum 
difference  in  temperature  between  the  upper  and  lower  sections  is  approximately  the  same.  Higher 
temperatures  are  known  to  stabilize  the  flame,  but  since  the  actual  temperature  at  extinction  varies  in 
the  experiment  in  a non-systematic  way,  the  uncertainty  in  LFL  and  CFR  caused  by  variations  in 
reactant  temperature  can  be  reduced  by  replication.  The  high  degree  of  repeatability  of  the  extinction 
conditions  indicate  that  such  random  errors  are  smaller  than  the  uncertainty  in  flow.  The  impact  of  the 
barometric  pressure,  which  ranges  between  about  98  kPa  and  100  kPa,  is  even  less  than  that  of  the 
initial  temperature  since  combustion  theory  predicts  almost  no  change  in  flame  stability  for  such  a small 
change  in  pressure. 

The  extinction  process  is  dynamic  and  the  response  time  of  the  MFCs  and  burner  are  non-zero. 
This  means  that  the  exact  conditions  at  the  nozzle  exit  and  in  the  flame  during  the  precise  point  of 
extinction  are  not  measured.  The  uncertainty  due  to  this  behavior  is  reduced  by  ensuring  that  changes 
in  flow  conditions  occur  at  a rate  slower  than  the  response  time  of  the  burner  and  control  system, 
which  is  about  10  s.  The  conditions  at  extinction  are  taken  to  be  the  readings  just  after  the  change  in 
flow  setting.  Thus,  the  uncertainty  can  be  taken  as  one  half  the  increment  between  the  previous  and 
final  step.  This  value  varies  among  tests,  but  is  usually  less  than  0.5  % of  the  recorded  stretch  rate  or 
mole  fraction. 

Intentionally  changing  the  burner  geometry  significantly  effected  the  measured  LFL,  as  was 
discussed  in  the  results  section.  This  is  distinct  from  small  changes  due  to  imperfections  in  the  burner 
or  misalignments  in  assembly.  The  burner  was  disassembled  for  cleaning  a number  of  times  and 
physically  relocated  from  one  laboratory  to  another.  No  extraordinary  care  was  taken  to  reassemble 
the  nozzles  precisely  in  the  same  manner  each  time.  It  is  estimated  that  the  nozzle  spacing  and 
centerline  alignment  could  have  varied  by  as  much  as  0.5  mm.  As  long  as  an  entire  test  sequence  was 
conducted  without  disassembling  the  burner,  no  additional  uncertainty  in  LFL  was  found  distinct  from 
the  random  errors  associated  with  run-to-run  variations. 
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The  variations  in  flammability  limit  (>  8 % of  the  LFL  for  methane,  more  for  R-32)  caused 
by  gross  changes  in  the  burner  geometry  (±  20  % in  spacing  and  diameter,  straight  nozzles  versus 
contoured  nozzles)  are  associated  with  non-idealities  in  the  flame  structure.  The  theory  is  premised  on 
the  following: 

• that  the  flame  is  axisymmetric; 

• that  the  radial  gradients  of  scalar  quantities  are  much  less  than  the  axial  gradients; 

• that  all  chemical  species  diffuse  at  equal  rates,  and  at  about  the  same  rate  as  the  diffusion  of 
heat  and  momentum;  and 

• that  the  system  is  adiabatic. 

The  first  two  assumptions  increase  in  validity  as  burner  diameter  and  spacing  become  large,  and  when 
the  initial  velocity  profile  is  flat.  The  diffusion  coefficients  vary  among  the  individual  species  by  an 
order  of  magnitude,  but  are  not  affected  directly  by  geometry.  Likewise,  the  gradients  in  the  primary 
reaction  zone  which  drive  diffusion  are  established  by  the  chemistry  and  not  the  exact  geometry. 
However,  the  rate  of  diffusion  of  all  species  in  the  radial  direction  is  lessened  as  the  flame  becomes 
more  one  dimensional. 

Heat  transfer  to  the  burner  walls  by  conduction  and  to  the  surrounding  environment  by 
radiation  causes  the  flame  to  be  non-adiabatic.  Conduction  losses  are  small,  in  general,  and  are 
lessened  as  the  distance  between  the  flame  and  the  burner  are  increased,  either  by  increasing  the  nozzle 
diameter  or  the  spacing. 

Radiation  losses  are  also  small  in  absolute  terms,  but  increase  directly  with  the  mean-beam- 
length  of  the  high  temperature  flame  zone.  Thus,  radiation  losses  are  enhanced  as  the  size  of  the 
burner  is  increased,  which  is  the  trend  opposite  from  the  first  three  assumptions  above  and  the 
conduction  losses.  It  is  possible  that  the  dichotomy  in  behavior  related  to  scale  is  responsible  for  the 
complicated  relationship  between  gross  changes  in  burner  geometry  and  the  measured  flammability 
limits. 

Fortunately,  the  LFL  of  methane  measured  in  the  current  study  with  nozzle  diameter  and 
spacing  near  12  mm  is  consistent  with  that  of  other  researchers  using  different  counter-flow  burner 
designs  in  normal  and  microgravity  experiments.  There  are  no  other  reported  measurements  of  the 
LFL  for  R-32  using  counter-flow  burners  of  any  design,  but  the  value  determined  in  the  current  12  mm 
burner  is  consistent  with  those  measured  using  the  ASTM  apparatus  by  numerous  people. 

Comparison  to  Chemical  Kinetics  Flame  Model 

The  extinction  of  dual,  one-dimensional,  counter-flow  hydrocarbon/air  flames  has  been  numerically 
modeled  with  some  success  by  a number  of  researchers  [14,15].  The  addition  of  fluorine  to  the 
hydrocarbon  kinetics  scheme  greatly  expands  the  number  of  molecular  species  that  need  to  be  tracked 
in  the  calculations.  The  computational  penalty  associated  with  including  a complete  chemical  kinetics 
mechanism  such  as  the  one  developed  by  Burgess,  et  al.  [7],  in  a two-dimensional  computation  that 
accounts  for  the  non-ideal  nozzle  flow  is  impractical,  and  may  be  unnecessary  to  explain  the  qualitative 
behavior  observed  in  the  experimental  methane  and  R-32  flames. 

As  a more  tractable  alternative,  the  flames  were  assumed  to  be  one-dimensional,  freely 
propagating,  and  adiabatic.  The  structure  of  such  a flame  can  be  calculated  in  a straightforward 
manner  using  the  PREMIX  code  [16]  developed  by  Sandia  National  Laboratories  with  the  chemical 
kinetics  package  CFIEMKIN  [17].  The  methane/air  chemistry  was  based  upon  the  GRI  mechanism 
[18],  while  the  F/C/H/O  mechanism  developed  by  Burgess  et  al.  [7]  was  used  to  model  the  detailed 
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fluorine  chemistry  with  reactions  up  through  C2.  A total  of  780  chemical  reactions  were  assumed  to 
take  place  among  83  species,  as  listed  in  Appendix  D. 

The  PREMIX  flame  code  has  a number  of  options  which  control  the  calculation  procedure  and 
can  impact  the  predicted  species  mole  fractions,  temperature  profile,  and  flame  speed.  In  the  current 
study  the  secondary  effect  of  the  temperature  gradient  on  mass  diffusion  (Soret  effect)  was  included. 
Upwind  differencing  was  used  for  the  convective  term,  and  the  flame  was  assumed  to  be  anchored  at 
the  location  where  the  temperature  reached  127  °C.  The  parameters  GRAD  and  CURV  control  the 
development  of  the  grid  spacing,  with  small  values  of  each  restraining  the  maximum  first  and  second 
derivatives  in  the  species  profiles  that  will  be  tolerated.  The  former  was  set  at  0.1  and  the  latter  at  0.3. 
The  absolute  and  relative  tolerances  placed  on  convergence  of  the  Newton  iteration  were  l.OxlO'9  and 
l.OxlO"4,  respectively. 

The  required  number  of  grid  points  across  the  flame  in  a converged  solution  ranged  between 
126  for  a stoichiometric  CH4/air  flame  to  238  for  an  R-32/air  flame  approaching  its  lean  flammability 
limit.  The  computational  domain  extended  from  - 50  mm  to  + 1000  mm.  The  numerical  program  as 
received  from  Sandia  National  Laboratory  was  designed  to  run  on  a work-station,  but  the  code  was 
modified  for  the  current  effort  to  be  compatible  with  the  NIST  Convex  C3820  vector  machine. 
Depending  upon  the  initial  conditions  in  the  problem  and  the  accuracy  of  the  first  guess  for  the 
temperature  profile,  it  took  from  1000  s to  40,000  s of  CPU  time  to  reach  a converged  solution. 

Stoichiometric  methane  and  R-32  mixtures:  The  complete  combustion  of  a mixture  of  methane 
and  R-32  at  an  overall  equivalence  ratio  of  unity  is  given  by  the  following  expression: 

XCH2F2+  (1-x)  CH4+  (2-x)(02  +3.76  N2)  - C02+  2XHF+  2(l-x)H20  +3.76(2-X)N2  {R6} 

where  x is  the  mole  fraction  of  R-32  in  the  binary  fuel.  Computations  were  performed,  first,  in  the 
limiting  cases  of  x = 1 and  x = 0 to  compare  the  structure  of  a pure  R-32/air  flame  to  that  of  a 
methane/air  flame. 

The  upper  graph  in  Figure  19  is  a plot  of  the  temperature  in  the  two  flames.  The  initial 
temperature  at  the  left  is  25  °C.  In  the  methane  flame,  the  temperature  rises  steeply  within  the  first 
millimeter  and  reaches  1960  °C  by  the  end  of  the  computational  domain.  The  temperature  in  the  R-32 
flame  builds  up  more  slowly  but  eventually  reaches  almost  the  same  value  (1930  °C).  By  transforming 
the  distance  scale  to  a time  scale,  the  temperature-time  gradient  can  be  used  to  accentuate  the  difference 
in  temperature  build  up  within  each  of  the  two  stoichiometric  flames,  as  seen  in  Figure  20.  The 
magnitude  of  the  methane/air  peak  is  25-fold  greater  than  the  magnitude  of  the  R-32/air  peak,  and  the 
methane  peak  occurs  almost  an  order-of-magnitude  earlier  in  time. 

The  lower  portion  of  the  graph  in  Figure  19  compares  the  mole  fractions  of  CO,  OH  and  CH2  (ground- 
state)  in  the  CH4  and  CH2F2  flames.  The  carbon  monoxide  begins  to  form  earlier  in  the  R-32  flame, 
but  the  rate  of  formation  of  CO  in  the  CH4  flame  accelerates  and  reaches  its  maximum  mole  fraction 
sooner,  followed  by  a decay  to  the  final  equilibrium  value.  The  peak  CO  mole  fraction  is  about  the 
same  in  both  flames.  The  OH  mole  fraction  is  indicative  of  the  size  of  the  chain-propagating  radical 
pool,  and  is  shown  to  peak  in  the  R-32  flame  beyond  the  CO.  The  level  of  OH  is  about  an  order-of- 
magnitude  smaller  than  the  OH  in  the  CH4  flame.  The  ground-state  triplet  methylene  (CH2)  behaves  in 
a way  representative  of  other  small  hydrocarbon  radicals  (e.g.,  CH3,  CH).  It  reaches  a peak  at  a 
location  close  to  that  of  the  CO,  and  then  practically  disappears  shortly  beyond  the  OH  maximum  for 
both  fuels.  The  mole  fraction  of  CH2  in  the  R-32  flame  is  over  ten  times  lower  than  in  the  methane 
flame. 

When  the  R-32  and  methane  are  combined  into  a single  flame,  the  calculated  normal  flame 
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MOLE  FRACTION  TEMPERATURE 


Figure  19.  Computed  temperature  and  species  profiles  through  one-dimensional  methane/air  and 
R-32/air  flames 


39 


TEMPERATURE-TIME  GRADIENT,  °C/ms 


Figure  20.  Computed  rate-of-temperature  increase  across  one-dimensional  flame,  comparing 
R-32/air  and  methane/air  mixtures 
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speed,  v0,  for  a stoichiometric  mixture  decreases  from  407  mm/s  to  67  mrn/s  as  the  mole  fraction  of 
refrigerant  is  increased  from  0 to  1.0.  The  calculated  flame  speeds  are  shown  in  Figure  21,  along  with 
the  final  flame  temperatures  and  maximum  OH  and  H mole  fractions.  The  flame  temperature  varies 
only  slightly  with  the  fraction  of  R-32,  whereas  the  OH  and  H drop  monotonically  to  less  than  15  % of 
their  initial  value  as  x is  increased  from  0 to  1 . 

The  calculated  flame  speed  for  pure  methane/air  is  in  agreement  with  measurements  by  Linteris 
and  Truett  [19],  the  lone  experimental  study  identified  in  which  the  speed  of  an  R-32/methane/air  flame 
has  been  explicitly  determined.  They  used  a premixed,  laminar  nozzle  burner  in  which  increasing 
amounts  of  R-32  were  added  to  an  initially  lean,  stoichiometric  or  rich  methane/air  flame.  The  flame 
speed  was  determined  from  Schlieren  photographs  of  the  flame  cone  angle.  The  solid  diamonds  plotted 
in  Figure  21  are  their  data.  The  measured  flame  speeds  at  x=  0 and  x=  0.18  are  within  the 
uncertainty  of  those  predicted  by  the  PREMIX  model.  (The  numerical  uncertainty  for  lean  flames  is 
estimated  to  be  ± 20  mm/s  based  upon  repeated  calculations  using  different  initial  temperature  profiles 
and  grid  control  parameters.)  For  the  highest  mole  fraction  of  R-32  studied  in  the  experimental  flame 
(x=  0.46),  the  overall  equivalence  ratio  was  about  1.2  (even  though  <I>  based  upon  the  methane/air 
ratio  was  0.9).  The  difference  in  equivalence  ratios  between  the  measurements  of  Linteris  and  Truett 
and  the  PREMIX  prediction  (in  which  3>=  1.0)  may  account  for  the  45  mm/s  discrepancy  in  flame 
speed.  This  was  not  confirmed  with  PREMIX  because  the  fluorine  mechanism  is  uncertain  in  rich 
mixtures. 

Lean  R-32/air  mixtures:  The  impact  of  equivalence  ratio  on  the  structure  of  the  pure  R-32/air 
flame  also  has  been  examined  numerically,  and  compared  to  the  impact  of  <£  on  the  methane/air  flame. 
Figure  22  is  a plot  of  the  final  temperatures  and  peak  OH  and  H mole  fractions  as  a function  of  4>.  As 
one  would  expect,  the  temperature  of  the  refrigerant  flame  decreases  continuously  with  decreasing  4>. 

It  is  noteworthy,  however,  that  the  final  temperature  in  the  R-32  flame  exceeds  that  of  the  methane 
flame  when  the  equivalence  ratio  is  leaner  than  0.90.  Of  great  significance  is  the  difference  in  behavior 
of  the  peak  OH  mole  fraction  for  the  two  fuels.  For  the  R-32/air  flame,  not  only  is  the  level  of  the  OH 
much  less,  but  also  the  shape  of  the  curve  is  qualitatively  different.  The  calculated  OH  mole  fraction  in 
the  methane/air  flame  drops  by  a factor  of  ten  as  changes  from  1.0  to  0.5,  and  H-atom  by  a factor  of 
100.  In  the  R-32/air  flame  the  OH  mole  fraction  actually  increases  slightly  as  the  flame  moves  from 
stoichiometric  to  <E>  = 0.8,  but  remains  close  to  0.001  over  the  entire  range  of  equivalence  ratios 
examined.  The  H-atom  mole  fraction  decreases  monotonically  with  O,  but  not  as  steeply  as  calculated 
for  the  methane/air  flame.  The  slower  change  in  OH  and  H-atom  mole  fractions  with  decreasing 
equivalence  ratio  contributes  to  the  ambiguity  in  the  defining  a precise  flammability  limit  for  R-32/air 
flames. 

Figure  22  compares  the  propagation  speeds  for  the  two  flames  at  different  values  of  O.  The 
normal  flame  speed  drops  slowly  with  decreasing  equivalence  ratio  in  the  R-32  flame,  eventually 
attaining  a value  of  36.7  mm/s  for  $ = 0.68.  The  numerical  model  computes  a value  of  29  mm/s  for 
an  unstretched,  adiabatic  methane/air  flame  at  an  equivalence  ratio  of  0.48. 

R-125/32:  The  chemistry  of  two-carbon  hydrofluorocarbon  oxidation  is  included  in  the  kinetics 
scheme  of  Burgess  et  al.  [7],  which  allowed  the  flame  speed  of  R-125/32  mixtures  to  be  estimated  in  the 
Phase  II  report  [5].  Figure  23  is  a reprint  of  those  results.  The  flame  temperature,  velocity,  and  OH 
concentration  for  the  stoichiometric  R-32/air  flame  are  used  to  normalize  the  parameters  as  the  mole 
fraction  of  R- 125  in  the  fuel  mixture  is  increased  from  0 to  0. 14.  The  overall  equivalence  ratio  is  kept  at 
1.0  based  upon  the  definition  in  Equation  (1)  (see  p.  27).  The  temperature  drops  only  slightly  as  the  R-125 
fraction  is  increased.  However,  the  OH  decreases  dramatically,  dropping  to  about  5 % of  the  pure  R-32 
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R-32  MOLE  FRACTION  IN  FUEL 


Figure  21.  Effect  of  replacing  methane  with  R-32  on  stoichiometric  flame  speed  (open  diamonds), 
equilibrium  temperature  (circles),  and  peak  mole  fractions  of  H (triangles)  and  OH 
(squares).  Experimental  flame  speeds  [17]:  filled  diamonds. 
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Figure  22 
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Computed  equilibrium  temperature  and  peak  OH  and  H radical  mole  fractions  as  a 
function  of  equivalence  ratio,  comparing  methane/air  to  R-32/air  flames. 
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Figure  23.  Effect  of  equivalence  ratio  on  computed  flame  speed  and  characteristic  time  for 

chemical  reaction,  compared  to  estimated  residence  time  in  experimental  flames  at 
extinction 
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Figure  24.  Calculated  flame  parameters  as  a function  of  the  amount  of  R-125  on  a stoichiometric 
R-125/32/dry  air  flame  initially  at  25  °C 
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flame  when  the  R-125  fraction  is  0.14.  The  flame  speed  at  this  R-125  level  is  less  than  40%  of  the  pure 
R-32  flame. 

Impact  of  relative  humidity  and  initial  temperature:  All  of  the  flame  calculations  presented  so  far 
have  taken  the  air  to  be  dry.  Figure  25  demonstrates  what  happens  when  moisture  is  added  to  the  air. 

The  saturation  pressure  of  water  at  25  °C  is  3.17  kPa.  This  means  that  when  the  air  pressure  is 
atmospheric  (101  kPa)  the  mole  fraction  of  water  at  100%  relative  humidity  is  0.031.  The  actual  mole 
fraction  of  water  used  in  the  calculations  is  reduced  by  a factor  equal  to  the  mole  fraction  of  air  in  the  dry 
air  plus  fuel  mixture.  The  upper  graph  in  Figure  25  is  for  a stoichiometric  10  % R-125/90  % R-32  flame, 
and  the  lower  graph  is  for  a lean  ($  = 0.70)  R-32/air  flame.  The  final  temperature  is  not  much  affected 
by  the  moisture,  but  the  peak  OH  mole  fraction  increases  almost  four-fold  in  the  R-125/32  flame. 
Interestingly,  the  normalized  flame  speed  is  enhanced  with  added  moisture  in  the  R-125/32  flame,  while  it 
is  reduced  in  the  R-32/air  flame.  This  seeming  inconsistency  can  be  explained  by  comparing  the  H-atom 
to  F-atom  ratio  in  the  two  flames.  The  pure  R-32/air  flame  has  an  H/F  ratio  of  1:1  independent  of  the 
stoichiometry.  This  means  that  there  is  no  excess  of  fluorine  atoms  to  tie  up  the  H-atoms  critical  to  the 
flame  propagation.  When  R-125  is  added  to  the  fuel,  the  H/F  ratio  drops  below  1:1  (0.826  for  the  10% 
R-125  flame).  Hence,  the  water  brings  the  OH  and  H levels  above  the  threshold  vital  to  maintaining  the 
combustion  reaction. 

The  final  parameter  investigated  was  the  initial  temperature.  Figure  26  (reprinted  from  the  Phase 
II  report  [5])  shows  how  the  velocity,  OH  level  and  final  temperature  are  impacted  as  the  initial  mixture 
temperature  is  increased.  The  system  is  a stoichiometric  10%  R-125/90  % R-32/dry  air  flame,  with  the 
values  at  25  °C  used  to  normalize  the  parameters.  An  initial  temperature  of  65  °C  has  an  imperceptible 
effect  on  the  final  normalized  temperature  and  a small  positive  effect  on  the  OH  mole  fraction.  The 
velocity  of  the  flame  increases  by  about  25  % . By  comparing  the  impact  of  the  temperature  change  to  the 
impact  of  changing  the  mole  fraction  of  R-125,  one  can  estimate  that  an  increase  in  initial  mixture 
temperature  of  10  °C  would  produce  an  absolute  increase  in  the  CFR  of  approximately  0.5  % for  this 
particular  flame  system. 

The  Importance  of  Flame  Speed  and  Damkohler  Number  as  Measures  of  Flammability 

The  counter-flow  burner  experimental  results  cannot  be  predicted  directly  from  the  PREMIX/ 
CHEMK1N  calculations  since  flame  stretch,  heat  loss  and  buoyancy  have  been  excluded.  Because 
these  natural  quenching  processes  are  absent,  the  numerical  code  predicts  a non-zero  flame  propagation 
rate  for  mixtures  leaner  than  the  experimental  flammability  limit.  Westbrook  [20]  suggested  that 
mixtures  with  one-dimensional,  adiabatic  flame  speeds  predicted  to  be  less  than  50  mm/s  are  beyond 
the  flammability  limit  from  a practical  standpoint.  Bui-Pham  et  al.  [21]  considered  a similar  criterion 
for  identifying  the  rich  flammability  limit  of  a methanol/CO/air  mixture,  and  found  it  to  correspond  to 
the  condition  where  the  rate  of  the  primary  chain  branching  reaction  ( H + 02  — > OH  +0)  is  equal  to 
the  rate  of  the  primary  chain  terminating  reaction  (H  + H02  -»  H2  4-  02),  a suggestion  originally  put 
forth  by  Law  and  Egolfopoulos  [22].  Using  50  mm/s  as  a qualitative  measure  of  the  flammability 
boundary,  then,  the  current  numerical  study  predicts  (see  Figure  23)  a practical  lean  flammability  limit 
of  0.53  for  CH4/air  and  0.77  for  R-32/air.  These  values  for  LFL  are  not  far  from  some  of  those 
measured  in  various  experimental  studies,  as  summarized  in  Tables  5 and  6. 

Chung  et  al.  [23]  demonstrated  that  extinction  is  likely  when  tf,  the  fluid  mechanical  residence 
time  in  the  flame,  is  less  than  the  characteristic  chemical  reaction  time,  tc.  The  Damkohler  number  is 
the  ratio  of  these  two  values,  D = xf  /xc,  so  that  a value  less  then  one  suggests  that  extinction  is  likely. 

A characteristic  fluid  residence  time  can  be  estimated  from  the  conditions  in  the  experimental  burner 
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100 


RELATIVE  HUMIDITY  IN  AIR  @ 25  °C 


Figure  25.  Effect  of  relative  humidity  on  (a)  10  % R-125/90  % R-32/air  flame  (0  = 1),  and 
(b)  R-32/air  flame  (0=0.7) 
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Figure  26.  Effect  of  initial  temperature  on  10  % R-125/90  % R-32/dry  air  flame  (4>  = 1) 
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Table  5.  Lean  flammability  limits  of  CH4/air  («25  °C)  using  different  experimental  methods 


Author(s) 

Method 

Conditions 

CH4/air  Lean 
Limit,  <J>0 

Richard  and  Shankland  [13] 

ASTM  E 681 

5 liter,  match 

0.51 

Zabetakis  [26] 

propagating  flame  tube 

extinction 

0.50 

this  work 

counterflow,  twin  flame 

linear  extrapolation 

0.49  ± 0.01 

Richard  and  Shankland  [13] 

ASTM  E 681 

5 liter,  match 

0.48 

Ishizuka  and  Law  [27] 

counterflow,  twin  flame 

linear  extrapolation 

0.48 

Yamaoka  and  Tsuji  [28] 

Tsuji  burner 

flame  location 

0.47 

Maruta  et  al.  [24] 

counterflow,  twin  flame,  pg 

turning  point 

0.47 

Sorenson  et  al  [29] 

coaxial  (tent)  flame 

flame  angle 

0.40 

Table  6.  Lean  flammability  limits  of  R-32/air  (~  25  °C)  using  different  experimental  methods 


Author(s) 

Method 

Conditions 

CH2F2/air 
Lean  Limit, 

$0 

Dekleva  et  al.  [12] 

5 cm  tube  (ICI) 

hot  wire 

1.11 

Richard  and  Shankland  [13] 

4 liter  tube 

match 

0.84 

Dekleva  et  al.  [12] 

ASTM  E 681,  5 liter 

hot  wire 

0.81 

Grob,  D.  [30] 

ASTM  E 681,  12  liter 

hot  wire 

0.81 

Richard  and  Shankland  [13] 

ASTM  E 681,  5 liter 

hot  wire 

0.79 

this  work 

counterflow,  twin  flame 

linear  extrapolation 

0.78  ± 0.04 

Dekleva  et  al.  [12] 

ASTM  E 681,  5 liter 

hot  wire 

0.77 

Dekleva  et  al.  [12] 

ASTM  E 681,  12  liter 

match 

0.77 

Dekleva  et  al.  [12] 

Autoclave,  8 liter 

hot  wire 

0.75 

Richard  and  Shankland  [13] 

ASTM  E 681,  5 liter 

spark 

0.74 

Ohnishi,  H.  [31] 

ASTM  E 681,  5 liter 

paper  match 

0.71 

Richard  and  Shankland  [13] 

ASTM  E 681,  5 liter 

match 

0.69 
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and  compared  to  the  characteristic  time  for  chemical  reaction  as  estimated  from  the  numerical 
simulation.  The  characteristic  fluid  residence  time  scales  with  the  distance  between  the  burner  outlet 
and  the  stagnation  plane,  divided  by  the  outlet  velocity  (i.e.,  the  inverse  of  the  global  stretch  rate, 

1/Kg).  The  residence  time  at  the  extinction  limit  as  measured  in  the  current  study  is  plotted  in  Figure 
23 . The  open  diamond  symbols  correspond  to  the  R-32/air  flame  and  the  filled  diamonds  refer  to  the 
methane/air  flame.  Also  plotted  in  Figure  23  (solid  triangles)  is  the  inverse  of  the  stretch  rate  near 
extinction  which  was  determined  by  Law  et  al.  in  their  counter-flow,  premixed  methane/air  burner, 
using  the  slope  of  the  local  velocity  on  the  center  line  in  the  preheat  zone  [3] . The  two  methane  flame 
data  sets  are  in  reasonable  agreement,  suggesting  that  the  global  stretch  rate  at  extinction  may 
reasonably  approximate  the  local  stretch  rate  at  extinction  for  the  lean  conditions  examined  in  the 
current  study,  a conclusion  also  reached  by  Maruta  et  al.  [24]  based  upon  the  work  of  Kobayshi  and 
Kitano  [25]. 

The  dotted  lines  shown  in  Figure  23  are  drawn  through  the  characteristic  reaction  times,  uc, 
determined  from  the  PREMIX/CHEMKIN  results  (symbolized  by  the  open  and  filled  squares  for  R-32 
and  methane  flames,  respectively).  The  reaction  time  is  approximated  by  the  transit  time  between  the 
location  of  the  flame  anchoring  temperature,  127  °C  (as  suggested  in  [17]),  and  the  position  of  the  peak 
H mole  fraction.  The  peak  in  H was  selected  as  a marker  for  reaction  time  because  of  the  importance 
of  H-atom  to  flame  propagation,  and  because  the  peak  was  found  at  a location  close  to  the  maximum 
levels  of  OH,  a radical  critical  to  the  burn-out  of  CO. 

By  comparing  the  numerically  calculated  tc  to  the  experimentally  determined  xf,  one  can  see 
that  extinction  in  the  actual  stretched  methane/air  flame  is  predicted  reasonably  well  by  the  PREMIX 
model  when  the  two  characteristic  times  are  about  equal  to  each  other;  i.e.,  D~  1.  The  same  cannot  be 
said  of  the  R-32/air  flame.  The  numerical  calculations  suggest  that,  for  all  equivalence  ratios  leaner 
than  0.9,  the  flame  should  be  more  robust  than  the  experimental  data  indicate.  This  discrepancy  may 
be  explained  three  possible  ways. 

A first  possibility  is  that  the  radiative  and  conductive  heat  losses,  which  are  not  included  in 
model,  are  more  significant  (and  therefore  more  detrimental)  to  the  refrigerant  flame  than  to  the 
hydrocarbon  flame.  While  radiation  heat  loss  has  been  included  by  others  in  PREMIX  models  of 
hydrocarbon  flames  [32],  an  estimate  of  the  relative  importance  of  the  radiation  in  the  two  flames 
studied  here  can  be  made  by  following  the  approach  of  Hertzberg  [33].  As  discussed  in  the  phase  II 
report  [5],  the  heat  loss  due  to  radiation  leads  to  a limiting  flame  speed  of  okp/rTf3/cpp,  where  o is  the 
Stefan  Boltzmann  constant,  kp  is  the  gray  gas  absorption  coefficient,  lr  is  the  radiation  length  scale,  Tf  is 
the  flame  temperature  in  Kelvin,  cp  is  the  specific  heat  of  the  flame,  and  p is  the  gas  density.  This 
limiting  flame  speed  can  be  compared  for  each  of  the  fuels  at  Kg  = 40  s'1,  corresponding  to  an 
extinction  equivalence  ratio  of  0.90  for  the  R-32/air  flame,  and  O = 0.52  for  methane.  At  these 
conditions,  the  calculated  flame  thickness  is  approximately  the  same,  but  the  flame  temperature  varies 
significantly:  2140  K for  R-32/air  and  1510  K for  methane/air.  While  there  is  some  difference  in  kp 
between  the  two  fuels,  it  about  cancels  with  the  change  in  density  due  to  temperature.  As  a result,  the 
radiation  flame  speeds  scale  with  Tf3,  whence  the  effect  of  radiant  heat  loss  is  almost  three  times  greater 
in  the  refrigerant  flame  than  in  the  methane  flame.  The  higher  temperature  in  the  R-32/air  flame  also 
leads  to  greater  heat  loss  to  the  cooled  burner  due  to  conduction,  which  scales  with  X AT//C  and  can  be 
significant  due  to  the  low  gas  velocities.  The  conduction  length  scale  lc  is  the  same  in  each  flame,  but 
the  product  of  the  thermal  conductivity,  X,  and  the  temperature  difference  between  the  flame  and  the 
burner,  AT,  is  just  nearly  twice  as  high  in  the  refrigerant  flame.  The  impact  of  heat  loss  on  extinction 
prediction  is,  thus,  more  significant  in  the  R-32 /air  flame,  and  adiabatic  calculations  of  flame  speed 
near  the  lean  limit  are  more  likely  to  over  predict  flame  speeds  when  R-32  is  the  fuel. 

A second  explanation  is  that  the  inverse  of  the  global  stretch  rate  is  a not  a good  indicator  of 
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the  fluid  residence  time  in  an  R-32/air  flame.  However,  it  is  not  necessary  for  the  absolute  value  of  the 
global  stretch  rate  to  be  precise,  only  that  it  vary  from  the  true  local  stretch  rate  by  a multiplicative 
constant.  The  validity  of  this  assumption  is  supported  by  the  agreement  between  the  LFL  measured 
with  the  counter-flow  burner  and  previously  reported  values,  and  by  the  work  done  in  [25].  Local 
velocity  and  temperature  measurements  through  the  flame  would  provide  a more  appropriate 
measurement  of  the  fluid  residence  time,  but  may  not  be  possible  under  conditions  approaching 
extinction. 

A final  consideration  is  that  the  chemical  kinetics  mechanism  is  incomplete  or  contains 
incorrect  rate  coefficients  for  the  fluorine-containing  reactions.  Linteris  and  Truett  [19]  found  the  same 
mechanism  adequate  to  predict  the  flame  speeds  in  their  premixed  R-32/methane/air  burner,  but  did  not 
attempt  to  model  the  system  with  a hydrogen/fluorine  ratio  less  than  3:1.  Considering  the  paucity  of 
flame  data  under  high  fluorine  loads  like  those  modeled  here,  a large  measure  of  uncertainty  remains  in 
the  chemical  kinetic  scheme. 

Non-linear  Extrapolation  to  the  Zero-Stretch  Condition 

The  basis  for  assuming  that  the  fundamental  LFL  can  be  obtained  by  a linear  extrapolation  of  the 
extinction  mole  fraction  to  a global  zero-stretch  condition  is  the  satisfactory  agreement  between  the 
experimental  measurements  and  a straight-line  fit  (e.g..  Figures  12  and  15).  An  identical  approach  was 
used  effectively  by  Wang  et  al.  [34]  for  pre-vaporized  benzene/air  mixtures.  Although  a strong 
correlation  is  undeniable  from  a statistical  analysis  over  stretch  rates  between  about  30  s'1  and  60  s'1, 
the  critical  lower  stretch  conditions  necessary  to  confirm  linearity  are  unattainable  in  the  burner  due  to 
the  dominance  of  buoyancy.  Recent  micro-gravity  experiments  have  shown  conclusively  that  the 
extinction  mole  fraction  of  methane  does  not  vary  in  a linear  fashion  as  the  stretch  rate  approaches  zero 
[24].  A more  detailed  analysis  of  the  reactions  is  necessary  to  better  understand  the  behavior  of  the 
flame  at  lower  stretch  rates. 

The  inverse  of  the  chemical  reaction  time,  ic,  is  a measure  of  the  rate  of  chemical  reaction  in 
the  burner.  At  extinction,  the  Damkohler  number  is  close  to  unity;  hence,  the  rate  of  reaction  is  of  the 
order  of  l/if,  which  is  to  say  Kg  at  extinction.  If  an  Arhenius  expression  is  assumed  for  the  rate 
coefficient,  and  the  reaction  rate  is  taken  to  be  first  order  in  fuel  and  oxygen  concentration,  then  the 
following  empirical  expression  can  be  used  to  model  the  extinction  stretch  rate: 

reaction  rate  oc  l/ic  « l/if  « Kg  = C,  [fuel]  [OJ  exp  (-C2/Tf),  (2) 

where  C,  is  a proportionality  constant,  C2  is  the  activation  temperature,  Tf  is  the  flame  temperature, 
and  the  brackets  indicate  concentration  in  moles  per  unit  volume. 

If  no  heat  losses  occur,  the  equilibrium  temperature  (in  Kelvin)  of  one  mole  of  R-32  in  air  can 
be  written  in  terms  of  the  adiabatic,  stoichiometric  equilibrium  temperature,  Tad,  the  initial  temperature, 
Tj,  and  the  equivalence  ratio,  assuming  the  excess  air  acts  as  a heat  sink  and  the  specific  heat  of  the 
mixture  per  unit  mass  is  unchanged: 

Tf  = Tj  + ( Tad-  Tj)  $ / [$+0,726  (1-$)]  (3) 

The  mole  fractions  of  R-32  and  oxygen  are  equal  to  $/($  + 4.76)  and  l/($+4.76),  respectively.  The 
concentration  of  the  R-32  and  oxygen  are  proportional  to  their  respective  mole  fractions  and  the  molar 
density,  which  decreases  with  increasing  temperature.  Thus, 
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[R-32]  [OJ  oc  <&/($  + 4.76  )2/  T2. 


(4) 


For  this  flame,  Tad  at  $ = 1.0  is  calculated  from  reference  [6]  to  be  2211  K.  The  empirical 
constants  are  found  by  matching  the  experiment  at  30  s'1  and  102  s'1.  The  relation  between  the 
extinction  stretch  rate  and  the  equivalence  ratio  becomes 

Kg($)  = 7.17  x 10H  [0/(0  + 4.76  )2  / T 2 ] exp  (-16  200/Tf)  , s'1,  (5) 

with 

Tf  = 373  + 1838  O / [O+0.726  (1-0)],  K.  (6) 

The  equivalence  ratio  can  be  written  in  terms  of  the  mole  fraction  of  R-32  in  the  mixture;  i.e.,  O = 

4.76  XR_32/(1-XR_32).  The  dash-dot  line  in  Figure  27  is  a plot  of  Equation  (5),  indicating  a far  from 
straight  region  for  stretch  rates  below  30  s'1.  The  circles  in  Figure  26  are  a re-plot  of  the  lean  portion  of 
the  R-32/air  experimental  data.  For  comparison,  the  solid  line  is  the  linear  fit  through  the  experimental 
data  with  stretch  rates  of  70  s'1  and  less.  The  linear  extrapolation  yields  a zero  stretch  mole  fraction  of 
0.131,  while  the  non-linear  extrapolation  predicts  a zero  mole  fraction  lower  limit  for  R-32  at  a zero 
stretch  condition. 

Equation  (5)  may  be  physically  based,  but  it  provides  no  useful  estimate  of  the  practical  lower 
limit  because  it  does  not  account  for  the  heat  loss  due  to  radiation.  Theoretical  analysis  of  the  flame 
structure  allows  one  to  identify  a radiation  limit  that  comes  into  play  at  very  low  stretch  rates  [35]. 

The  radiation  loss,  Q^,  in  Joules  per  kilogram  of  R-32,  from  a disk  shaped  flame  can  be  estimated  by 

Qrad=  SifOkpTf4  (m VP)  (20  + 4.76)/(  0.052$  ),  (7) 

where  lq,,  is  the  absorption  coefficient , a is  the  Stefan-Boltzmann  constant  and  P is  the  pressure. 
Equation  (7)  indicates  that  the  radiation  loss  increases  in  an  unbounded  fashion  directly  with  the  flow 
time;  hence,  as  Kg  (»  l/if ) approaches  zero,  the  chemical  reaction  will  be  quenched. 

Equation  (5)  still  holds  when  the  system  is  non-adiabatic,  but  the  flame  temperature  is  reduced 
by  the  absolute  value  of  the  ratio  of  Qrad  to  the  enthalpy  of  combustion  of  the  fuel,  or 

Tf,  rad  = T,  + (1  - | Qrad/AHJ  | ) {(  Tad,*=1  - T|)  $ / [$+0,726  (1-0)]}.  (8) 

The  absorption  coefficient  was  estimated  in  [5]  to  be  about  0.8  m'1,  AHC  = -9.35  MJ/kgR_32,  the  initial 
temperature  is  373  K,  and  P is  101  kPa.  The  flame  temperature  accounting  for  radiative  loss  is,  thus, 

Tf-  rad  = 373  + [l-  4.46xl0'16  (Tf>rad4/Kg  )(2o  + 4.76)/(0.052o)]{1838  O / [O+0.726  (1-0)]}.  (8') 

The  constants  in  Equation  (5)  can  be  recalculated  based  upon  the  radiation-corrected  temperature, 
yielding  C,  = 7.35xl013  K2 s1  and  C2  = 18  300  K.  Equations  (5)  and  (8')  can  be  solved  iteratively  to 
find  the  extinction  equivalence  ratio  (or  XR.32)  as  a function  of  Kg  for  a flame  with  radiative  heat  loss. 
The  dashed  line  in  Figure  27  shows  the  result.  Accounting  for  radiative  heat  loss  does  two  things;  first, 
it  shifts  the  extinction  mole  fraction  curve  upward,  and  second,  it  demonstrates  a true  lower  limit,  as 
indicated  by  the  * in  Figure  27.  No  solutions  are  mathematically  obtainable  for  XR.32<  0.118.  The 
non-adiabatic,  non-linear  theory  predicts  a lower  value  for  LFL  than  the  simple  linear  extrapolation; 
however,  other  affects  such  as  conduction  losses,  preferential  diffusion,  and  two-dimensional  flow  act 
to  move  the  practical  LFL  upward. 
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R-32  Mole  Fraction 


Figure  27.  Comparison  of  non-linear  and  linear  extrapolations  to  the  zero  stretch  condition  for 
R-32/air  mixtures  initially  at  100  °C 
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A similar,  non-linear  analysis  can  be  applied  to  the  Critical  Flammability  Ratio  calculation. 

For  simplicity,  we  have  examined  the  behavior  of  an  R-32/N2  mixture  first.  Assume  that  the  air/R-32 
ratio  is  constant  and  equal  to  the  most  robust  condition  (<5  = 1.18)  as  an  increasing  amount  of  nitrogen  is 
added.  The  mole  fraction  of  the  nitrogen  in  the  R-32  plus  N2  mixture  is  x-  The  concentrations  of  the  R- 
32  and  oxygen  are  diluted  by  the  N2  such  that  they  can  be  expressed  as 

[R-32]  [OJ  oc  [(1-  x)/(5.03-4.03X)]2  / Tf2  (9) 

The  nitrogen  acts  as  a heat  sink  to  reduce  the  adiabatic  flame  temperature  according  to  the  relation 

Tf  = Tj  + ( Tad-T;)  [6.19(1-  X)/(6.19-5.19x)].  (10) 

The  adiabatic  temperature  increase,  ( Tad-  Tj),  is  reduced  from  1838  K to  1732  K based  upon  the 
equilibrium  value  computed  with  reference  [6]  for  <5  = 1.18  when  x=0.  The  radiative  loss  can  be 
calculated  if  it  is  assumed  that  the  absorption  coefficient  remains  about  constant.  Then, 

Qrad=3rfokpTf4  (RTj/P)  (6. 19-5. 16X)/[0. 052(1 -X)];  (11) 

and,  since  the  N2  does  not  affect  the  enthalpy  of  combustion, 

Tf,  rad  = 373  + 

{1-  4.46X10-16  (Tf  ) (6.19-5. 19X)/[0.052(1-X)]}{1732[6.19(1-  x)/(6.19-5.19x)]},  K.  (12) 

Figure  28  is  a plot  of  the  extinction  stretch  rate  measured  with  nitrogen  added  to  the  R-32.  The 
open  circles  are  the  experimental  data,  taken  with  the  temperature  equal  to  100  °C  and  a relative 
humidity  of  50  % (referenced  to  23  °C).  The  measured  extinction  values  for  x drop  for  stretch  rates 
below  30  s'1,  a phenomena  also  observed  with  R-125/R-32  mixtures  (Figure  16).  This  is  due  to 
buoyant  distortions  in  the  flame  at  low  stretch  rates.  The  values  for  C{  and  C2  for  R-32/N2  mixtures 
can  be  recalculated  from  the  measured  extinction  stretch  rate  for  x equal  zero  and  for  x close  to  the 
maximum  attainable.  The  extinction  stretch  rate  becomes 

Kg(x)  = 5.32xl014{[(l-x)/(5.03-4.03X)]2  / T } exp  (-21  600/Tf>rad),  s'1.  (13) 

Equation  (13)  is  plotted  as  the  dashed  line  in  Figure  28.  The  equation  predicts  extinction  for  x ~ 0.64 
(indicated  by  the  asterisk),  which  is  a little  greater  than  the  value  of  0.603  obtained  by  a linear  fit 
(dotted  line)  of  all  the  experimental  data. 

An  attempt  was  made  to  model  the  inhibiting  effect  of  the  R-125  as  with  nitrogen  (i.e., 
assuming  extinction  is  the  result  of  dilution  and  heat  absorption).  The  similarity  in  shape  of  the  x vs. 

Kg  curves  for  N2  and  R-125  argues  for  this  approach.  However,  the  first-order  relation  between  the 
fuel  (i.e.,  R-32)  and  reaction  rate  did  not  yield  the  close-to-straight-line  fit  observed  in  the  experiments 
for  meaningful  values  of  activation  temperature.  In  an  alternative  approach  the  R-125  was  treated  as  a 
fuel  that  contributes  to  the  heat  release  and  maintains  the  overall  stoichiometry  constant  as  x is 
increased.  The  activation  temperature  would  be  expected  to  be  higher  due  to  the  much  lower  reactivity 
of  the  R-125  when  compared  to  R-32.  However,  the  simple  model  used  to  predict  the  R-32  LFL  and 
the  CFR  of  nitrogen/R-32  assumes  that  a single  value  of  C2  exists  over  the  entire  range  of  x-  This  is 
invalid  for  a chemically  acting  substance,  and  may  be  the  reason  that  the  second  attempt  to  fit  the 
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Figure  28.  Comparison  of  non-linear  and  linear  extrapolations  to  the  zero  stretch  condition  for 
R-32/N2/air  mixtures  initially  at  100  °C 
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Figure  29.  Adiabatic,  equilibrium  temperature  and  mole  fractions  of  H,  H2,  and  OH  for 
R-245ca/R-32  mixtures  with  dry  air  initially  at  298  K. 
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Adiabatic  Equilibrium  Temperature 


experimental  data  for  R-125/32  mixtures  to  the  form  of  Equation  (2)  was  unsuccessful. 

In  the  case  of  R-245ca,  previous  unpublished  measurements  by  others,  using  the  ASTM 
constant  volume  apparatus,  identified  a region  of  flammability  in  humid  air  ignited  at  100  °C,  but  none 
when  the  initial  temperature  was  close  to  the  ambient.  Temperature  sensitive  and  somewhat  ambiguous 
results  were  also  found  in  the  current  research  program.  No  flame  could  be  sustained  for  pure  R- 
245ca/air  mixtures  in  the  counter-flow  burner  even  when  the  reactants  were  humidified  and  preheated 
to  100  °C.  A CFR  just  under  unity  was  identified  for  R-245ca/32  mixtures  by  extrapolating  the 
extinction  mole  fraction  to  zero-stretch  conditions  when  the  reactants  were  about  50  °C,  while  a small 
but  positive  extinction  stretch  rate  was  identified  as  being  attainable  for  pure  R-245ca/humid  air 
mixtures  preheated  to  100  °C. 

As  with  the  R-125/32  mixtures,  an  attempt  to  model  the  influence  of  Xr-245  on  the  extinction 
stretch  rate  was  unsuccessful  for  R-245ca/32  mixtures  because  of  the  complex  flame  chemistry.  In  this 
case,  the  heavier  refrigerant  inhibits  the  flame  in  low  concentrations  but  also  contributes  to  the 
flammability  when  it  is  the  dominant  fuel. 

The  adiabatic  equilibrium  temperatures  and  the  concentrations  of  H-atom,  H2  and  OH  in  R- 
245ca/32  mixtures  with  an  overall  equivalence  ratio  of  1.25  were  estimated  using  reference  [6],  and  are 
plotted  in  Figure  29.  A small  decrease  in  temperature  is  noted  when  small  amounts  of  R-245ca  are  first 
added  to  the  mixture,  but  the  equilibrium  temperature  is  almost  constant  for  Xr-245  between  0.1  and  1.0. 
Hence,  the  inhibiting  qualities  of  the  R-245ca  are  less  likely  due  to  heat  absorption,  and  more  likely 
associated  with  the  rapid  drop  in  key  intermediate  species  also  observed  in  Figure  29.  Although  the 
concentrations  in  the  flame  are  not  at  their  equilibrium  values,  the  high  diffusivity  of  atomic  and 
molecular  hydrogen  tends  to  smooth  out  the  gradients,  so  that  significant  changes  in  equilibrium  values 
of  H and  H2  are  reflected  in  similar  changes  within  the  flame.  The  relevance  of  these  equilibrium 
calculations  is  that  they  bolster  our  confidence  that  a linear  extrapolation  of  the  extinction  stretch  rate  to 
mole  fractions  of  R-245ca  below  those  that  can  be  obtained  in  the  counter-flow  burner  is  meaningful. 


CONCLUSIONS  AND  RECOMMENDATIONS 

This  research  has  demonstrated  that  a counter-flow  burner  is  well  suited  for  revealing  the  structure  of 
refrigerant/air  flames.  The  flame  is  stable  enough  to  investigate  the  relationship  between  equivalence 
ratio  and  the  stretch  rates  necessary  for  extinguishment.  A plot  of  the  mole  fraction  of  refrigerant 
versus  the  global  stretch  rate  shows  the  behavior  to  be  close  to  linear  down  to  velocities  where 
buoyancy  begins  to  distort  the  flame  and  heat  losses  become  significant. 

Summary  of  LFL  and  CFR  Measurements 

The  following  values  of  LFL  and  CFR,  with  their  uncertainties,  are  recommended  for  the  refrigerants 
examined  using  the  premixed,  counter-flow  research  burner  (12  mm  ± 0.2  mm  diameter  contoured 
nozzles,  spaced  12  mm  ± .0.2  mm  apart): 

R-32/air  LFL  (humidity  ratio  < 10'2  mg^yg^^,  T=  30  °C  ± 5 °C,  P = 98  kPa  ± 2 kPa): 

Xr32=  0.14  ± 0.006  (Note:  Uncertainty  reported  in  DOE/CE/238 10-98  improperly  as  ± 0.004) 

R-32/air  ETFL100  (humidity  ratio  = 8.5  mg^^/g^, ^ ± 0.5  mgwaler/gdr>. air , T = 100  °C  ± 10  °C,  P=  98 
kPa  ± 2 kPa): 

XR32=  0.13  ± 0.004 
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R-125/32  CFR  (humidity  ratio  < 10'2  mg^/g^ air,  T=  30  °C  ± 5 °C,  P = 98  kPa  ± 2 kPa): 

Xri25  = 0.18  ± 0.004 

R-125/32  CFR@100  (humidity  ratio  < 8.5  mg^/g^ * ± 0.5  mg^/g^ air,  T=  100  °C  ± 5 °C,  P = 
98  kPa  ± 2 kPa): 

Xri25  = 0.22  ± 0.004 

R-245ca/air  ETFL100  (humidity  ratio  = 8.5  mg^Jg^^  ± 0.5  mg T = 100  °C  ± 5 °C,  P=  98 
kPa  ± 2 kPa): 

Xr245=  1.35  ± 0.15 

R-245ca/air  LFL  (humidity  ratio  = 8.5  mg^/g^ * ± 0.5  mg^/g^^,  T=  50  °C  ± 5 °C,  P = 98 
kPa  ± 2 kPa): 

Non-flammable 


R-245ca/32  CFR  (humidity  ratio  = 8.5  mgwater/gdry air  ± 0.5  mg^Jg^ air,  T=  50  °C  ± 5 °C,  P = 98 
kPa  ± 2 kPa): 

Xr245  = 0.98  ± 0.02 


R-245ca/air  LFL  (humidity  ratio  < 10'2  mg^^/g^ air,  T=  100  °C  ± 5 °C,  P = 98  kPa  ± 2 kPa): 
Non-flammable 

R-245ca/32  CFR  (humidity  ratio < 10'2  v^g^^Jg^^,  T=  100  °C  ± 5 °C,  P = 98  kPa  ± 2 kPa): 
XR245  < 0.98 


N2/R-32  CFR@100  (humidity  ratio  = 8.5  mgwater/gdI>  ^ ±0.5  mg^/g^^,  T=  100°C±5°C,P  = 
98  kPa  ± 2 kPa): 

Xn2  = 0.61  ± 0.004 


These  results  demonstrate  that  the  major  objective  of  the  research  has  been  met:  to  design  and  operate 
a counter-flow  burner  that  produces  accurate  and  repeatable  flame  limits  and  critical  flammability  ratios 
at  temperatures  up  to  100  °C,  with  flammable  and  marginally  flammable  liquid  and  gaseous 
refrigerants. 


Industrial  Burner  Design 

The  measured  LFL  is  sensitive  to  some  geometric  variations,  the  most  significant  being  that  smaller 
diameter  burners  tend  to  widen  the  flammability  limits,  and  larger  diameter  burners  decrease  them. 

The  12  mm  diameter  burner,  however,  gives  results  consistent  with  LFL  measurements  of  R-32/air, 
methane/air  and  propane/air  mixtures  in  fixed  volume  experiments  and  flame  tubes.  It  is 
recommended  that  the  inner  diameter  of  the  counter-flow  burner  be  held  between  10  mm  and  12  mm, 
and  that  the  spacing  be  maintained  about  equal  to  the  diameter,  ± 1 mm.  The  exact  character  of  the 
nozzle,  the  nitrogen  annulus,  and  precise  alignment  of  the  burner  were  found  to  be  of  lesser  importance 
in  determining  the  zero-stretch  mixtures  as  long  as  no  changes  in  geometry  or  operational  procedure 
occur  within  a single  test  sequence. 

A new  burner  has  been  designed  based  upon  the  above  findings.  The  objective  has  been  to 


58 


MGKAKT  & AIR 


ippe  amxT 
WET 


TUBE 


24.09 

[612] 


IFPER  QUANT 

rnunmi 


/ 


ipts  aruNT 

OLfTLET 


neoaiPif 


DUTIET 


UO  COOLANT 
MJET 


Figure  30.  Overall  view  of  industrial-grade  counter-flow  burner 
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minimize  the  investment  in  facility  development  and  operation  to  the  refrigeration  industry,  while 
maintaining  a satisfactory  level  of  confidence  and  reproducibility  of  the  data.  Figure  30  shows  an 
overall  view  of  the  new  burner.  The  upper  and  lower  burner  tubes  are  similar,  and  made  of  three 
concentric  stainless  steel  tubes  of  standard  diameters  (1/2",  3/4"  and  1").  The  refrigerant  and  air  pass 
through  the  inner  tube;  the  19  mm  (3/4")  tube  directs  water  to  the  tip  of  the  burner  to  maintain  close 
control  of  the  gas  temperature  and  to  protect  the  outer  tube  from  overheating  in  the  hot  exhaust  stream. 
Four  adjustable  spacer  rods  permit  the  burner  tubes  to  be  aligned  conveniently.  A 152  mm  diameter 
Plexiglas  tube  shields  the  flame  and  directs  the  toxic  combustion  products  to  the  exhaust.  To  simplify 
the  design,  no  provision  is  made  for  a concentric  flow  of  nitrogen  gas  at  the  burner  outlets.  A heated 
electrical  wire  or  miniature  torch  can  be  rotated  onto  the  centerline  to  ignite  the  refrigerant/air  mixture. 
Thermocouples  are  located  in  the  water  at  the  burner  tip  to  monitor  the  initial  gas  temperature. 

Appendix  E contains  detailed  drawings  of  the  burner  components  and  a suggested  arrangement 
for  metering  all  the  flows.  Estimates  for  fabricating  a burner  ranged  between  $2500  and  $10,000  for 
the  first  unit.  It  is  recommended  that  several  be  built  and  delivered  to  different  laboratories.  A round- 
robin  testing  program  is  required  to  determine  if  the  industrial  burner  is  suitable  as  designed  to  produce 
consistent  results.  If  so,  and  if  the  operation  is  straight-forward,  consideration  should  be  given  to 
supplementing  or  replacing  the  constant  volume  ASTM  E 681-1994  method  with  a procedure  based 
upon  the  counter-flow  burner. 

Unresolved  Issues  and  Recommendations  for  Further  Research 

The  counter-flow  burner  test  method  circumvents  the  most  troubling  issue  that  is  intrinsic  to  ASTM  E 
681-1994:  Does  that  which  is  observed  in  the  ASTM  test  constitute  a self-propagating  flame?  The 
answer  to  the  question  is  tied  inextricably  to  the  details  of  the  design  of  the  ASTM  test  set-up,  the 
operational  procedure,  the  opinion  of  the  observer,  and  the  purpose  of  the  test.  The  ASTM  method 
combines  a fundamental  property  measurement  (which  should  be  unassailable)  with  an  assessment  of 
risk  (which  inherently  involves  personal  judgment).  Because  the  ignition  process  and  flame 
propagation  are  transient  and  three-dimensional  in  the  ASTM  apparatus,  meaningful  analysis  and 
extrapolation  to  other  environments  becomes  impossible. 

The  counter-flow  burner  method  is  based  upon  flame  extinction  (an  unequivocal  event)  to 
define  the  flammability  limit.  It  is  designed  to  measure  a fundamental  property,  divorced  from  the 
question  of  risk.  The  quasi-steady,  one-dimensional  nature  of  the  containerless  flame  makes  it 
amenable  to  detailed,  theoretical  analysis,  which  is  essential  before  risks  can  be  predicted  for  different 
possible  fire  scenarios. 

To  be  most  useful  to  the  refrigeration  industry,  a methodology  needs  to  be  developed  to  assess 
the  risk  of  a fire  from  a given  working  fluid  in  a particular  realistic  application.  If  one  knows  the 
laminar  flame  speed  of  a refrigerant/air  mixture,  and  how  it  varies  with  environmental  conditions,  then 
meaningful  predictions  of  the  fire  risk  associated  with  a given  geometry,  ignition  source,  and  leak 
scenario  can  be  made.  A number  of  carefully  planned  full  scale  tests  with  refrigeration  machine  mock- 
ups  could  be  performed  using  a range  of  class  2 and  class  3 refrigerants  to  bracket  the  uncertainty  of 
the  predictive  method. 

A comprehensive  model  of  the  counter-flow  burner  is  currently  lacking.  Further  research  is 
necessary  to  produce  such  a model,  which  would  then  permit  the  laminar  flame  speed  to  be  determined. 
The  following  tasks  are  recommended: 

• extend  the  chemical  kinetics  data  set  to  include  additional  refrigerants  of  interest,  such  as  ammonia 
and  the  three  carbon  propane  analogues; 
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• perform  a parametric  study  to  identify  the  primary  reaction  pathways  that  govern  the  flame 
extinction  process; 

• develop  a two-dimensional  fluid  mechanical  model  of  the  counter-flow  burner  that  accounts  for  heat 
loss  and  buoyancy; 

• couple  the  fluid  mechanical  model  to  the  chemical  kinetics  mechanism  and  conduct  a laboratory 
study  on  the  counter-flow  burner  to  demonstrate  the  model's  ability  to  predict  flame  extinction  over 
a range  of  burner  diameter,  nozzle  spacing,  and  inlet  conditions; 

• using  the  predictive  model  to  identify  worst  case  conditions,  design  a full-scale  test  to  determine  the 
fire  hazard  posed  by  a range  of  class  2 and  class  1/2  refrigerant  leaks. 

Taken  together  with  the  round-robin  testing  on  the  industrial  burner  recommended  above,  the 
results  of  this  research  would  allow  the  air-conditioning  and  refrigeration  industry  to  move  away  from 
an  ad  hoc  and  confusing  prescriptive  test  method  towards  a scientifically  sound,  performance-based 
safety 
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APPENDIX  A:  Step-by-step  Operation 


Special  Safety  Precautions: 

Make  sure  acid  resistant  gloves  are  worn  when  lighting  the  burner,  handling  the  scrubber  collection 
vats,  and  at  all  times  when  reaching  inside  the  chemical  hood. 

If  gas  flows  have  been  left  running  for  more  than  a few  minutes  when  the  flames  were  not  lit,  or  if 
there  is  any  reason  to  believe  an  excess  of  fuel  may  have  built  up  inside  the  burner  tube,  DO  NOT 
attempt  to  ignite  the  burner.  Instead,  shut  off  gas  flows  and  turn  on  the  blower  for  a few  minutes 
before  proceeding. 


Operation  Procedure  for  R-32,  R-125/32,  and  R-245ca/32  Mixtures: 

(Skip  steps  1-3  if  not  testing  R-245ca) 

1 . Set  up  the  liquid  delivery  system.  Place  the  syringes  into  position  on  the  pump,  inside  the  cooling 
coils,  and  attach  a thermocouple  to  the  end  of  each  syringe.  Connect  the  heated  delivery  lines  to  the 
burner  inlets  on  the  top  and  bottom  of  the  burner.  Turn  on  the  recirculating  cooling  bath  and  set  the 
temperature  to  approximately  0 ° C.  Turn  on  the  variacs  controlling  the  heating  tape  for  the  boilers  and 
the  heated  delivery  lines.  Set  the  temperature  of  the  boiler  at  least  30  ° C above  the  boiling  point  of  the 
liquid.  Set  the  delivery  lines  to  at  least  20  ° C above  the  boiling  point  of  the  liquid. 

2.  Monitor  the  temperature  of  the  syringes.  When  the  temperature  is  sufficiently  low  enough  to 
maintain  the  liquid  status  of  the  fuel,  the  syringes  may  be  filled.  For  R-245ca,  the  syringes  were  kept 
below  8.5  0 C to  maintain  a constant  density.  To  fill  the  syringes,  first  switch  the  valves  to  the  fill 
position.  Next,  push  the  plungers  all  the  way  forward  to  expel  any  air  in  the  syringes.  Then  invert  the 
bottle  of  liquid  fuel  and  attach  its  outlet  to  one  of  the  syringe  filling  tubes.  Open  the  valve  on  the  bottle, 
allowing  the  liquid  to  flow  into  the  syringe.  The  plunger  will  slowly  be  pushed  back  by  the  liquid. 

When  the  syringe  has  been  filled  to  the  desired  amount,  close  the  valve  on  the  bottle.  Then  turn  the 
syringe  valves  to  the  off  position.  Disconnect  the  bottle  and  repeat  procedure  for  the  other  syringe. 

3.  Turn  on  computer  and  open  flow  controller  program,  LIQUID.VI. 

(Skip  step  4 if  testing  R-245ca) 

4.  Turn  on  computer  and  open  flow  controller  program,  TWINFLAM. 

5.  Turn  on  main  water  valve  and  set  water  flow  rate  to  about  150  ml/min  through  the  heat  exchanger. 
Next,  turn  on  the  variacs  that  control  the  burner  heating  tape  and  the  heat  exchanger.  Set  the  heat 
exchanger  variac  to  an  initial  power  setting  of  85,  and  the  heat  tape  variac  to  an  initial  power  setting  of 
65.  This  will  sufficiently  pre-heat  the  burner  for  a 100  °C  testing  condition. 

6.  Check  the  gas  lines  running  to  and  from  the  bubblers,  to  make  sure  bubblers  are  connected.  Check 
the  that  the  bubbler  bath  is  cooled  to  13°C  for  50%  relative  humidity. 
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7.  Open  the  manual  valves  located  after  the  mass  flow  controllers  (MFCs)  and  open  all  gas  bottles.  At 
the  digital  control  box,  set  all  mass  flow  controllers  to  “ON”. 

8.  For  each  flow  controller,  set  the  digital  readout  at  the  control  box  to  the  initial  setting  at  which  it  was 
calibrated  (different  for  each  flow  controller).  With  the  gases  flowing,  adjust  the  gage  pressure  on  each 
bottle  to  140  kPa.  This  ensures  the  calibration  corrections  will  be  accurate. 

9.  Turn  off  the  gas  flow  of  the  fuel  at  the  digital  control  box.  While  the  air  is  still  flowing,  measure  the 
relative  humidity  and  dew  point  with  the  hand-held  hygrometer.  Record  these  values  for  later  use. 

10.  Determine  a flame  ignition  condition  (equivalence  ratio  and  stretch  rate)  for  the  pure  gas  fuel  and  air 
mixture  (without  R-245ca).  This  can  be  acquired  using  knowledge  of  the  lighting  conditions  from 
previous  tests  with  the  gas  fuel.  This  allows  for  easy  lighting,  then  the  vaporized  liquid  is  added  later. 

1 1.  In  the  flow  controller  program  (TWINFLAM  or  LIQUID. VI)  select  a file  name  under  which  the  data 
will  be  recorded,  and  set  the  test  number.  In  the  fuel  selection  box,  choose  the  liquid  fuel  and  gas  fuel  to 
be  used.  Type  in  the  measured  humidity,  and  nitrogen  co-flow  desired.  Also  type  in  the  initial  strain  rate 
and  equivalence  ratio  settings  for  the  lighting  conditions.  Set  the  gas  fuel  amount  to  100%  for  the  initial 
lighting  phase,  and  set  the  flow  controller  buttons  to  “ON”.  Once  all  selections  have  been  made,  start  the 
gas  flows  by  clicking  on  the  “Go”  arrow  in  the  upper  left  comer  of  the  screen. 

12.  Using  the  continuous  flame  igniter,  light  the  twin  flames.  Be  sure  to  wear  safety  gloves  during  this 
step,  and  whenever  reaching  inside  the  hood  since  HF  acid  may  be  present.  When  the  flames  are  lit,  turn 
on  the  blower  and  the  water  flow  for  the  scrubber.  Set  the  water  pressure  to  69  kPa.  Make  sure  the 
scrubber  collection  container  is  in  place. 

13.  The  flow  program  calculates  the  correct  setting  for  the  nitrogen  co-flow  rotameters.  Turn  on  the 
nitrogen  bottle  and  set  rotameters  as  indicated  by  the  program.  This  must  be  done  after  flames  are  lit. 
Flames  will  not  ignite  with  nitrogen  present. 

(Skip  steps  14-16  if  not  testing  R-245ca) 

14.  Once  flames  have  stabilized,  the  vaporized  liquid  fuel  may  be  added.  Set  the  pump  to  the  desired 
flow  rate,  and  enter  that  liquid  flow  rate  into  the  controller  program.  Turn  the  syringe  valves  to  the 
burner  delivery  position  and  start  the  pump.  Click  on  the  “Read  Data”  button  on  the  screen,  then  click  on 
the  “Go”  arrow.  The  program  will  calculate  the  vapor  flow  rate,  volume  percent  of  the  fuel  mixture, 
effective  equivalence  ratio  and  the  resulting  strain  rate. 

15.  Occasionally,  the  addition  of  the  liquid  fuel  will  cause  either  the  stretch  rate  or  effective  equivalence 
ratio  to  become  too  high,  resulting  in  a mixture  that  is  beyond  its  upper  flammability  range.  If  this 
happens,  repeat  steps  10  through  14  using  different  settings  for  the  lighting  condition  and/or  the  liquid 
flow  rate  until  a flame  can  be  sustained. 

16.  The  flame  is  now  a mixture  of  both  the  gas  fuel  and  the  vaporized  liquid  fuel.  The  percent  of  gas 
fuel  is  then  slowly  reduced  in  an  attempt  to  bum  the  vaporized  liquid  fuel  by  itself.  The  composition  of 
fuel  mixture,  and  either  the  strain  rate  or  equivalence  ratio  may  be  controlled.  The  other  parameter  will 
vary  accordingly.  Start  with  the  condition  resulting  from  the  addition  of  the  liquid  vapor.  Click  the 
“Read  Data”  button  to  record  the  initial  flows.  Select  whether  to  control  the  stretch  rate  or  the 
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equivalence  ratio  by  flipping  the  switch  to  point  towards  one  or  the  other.  Then  reduce  the  amount  of 
gas  fuel  by  setting  the  percent  gas  fuel  to  slightly  below  the  current  amount.  For  best  results,  make  very 
small  changes  to  the  flame  settings,  since  large  changes  to  the  gas  composition  may  cause  the  flame  to 
become  unstable  and  extinguish  prematurely.  Typically,  the  gas  fuel  is  reduced  in  steps  of  1%  of  the 
total  fuel  flow,  or  less.  Once  the  new  settings  have  been  selected,  activate  the  gas  flows  by  clicking  on 
the  “Go”  arrow  in  the  upper  left  comer.  Either  the  stretch  rate  or  equivalence  ratio  (whichever  has  been 
chosen  to  be  controlled)  is  held  constant,  while  the  other  adjusts  to  the  new  fuel  composition. 

(Skip  step  1 7 if  testing  R-245ca) 

1 7.  After  the  flame  stabilizes  at  the  chosen  condition,  new  settings  can  then  be  selected.  The  gas  fuel 
may  continue  to  be  decreased  until  it  is  eliminated.  Additionally,  the  stretch  rate  or  equivalence  ratio 
may  be  changed  from  its  initial  condition.  First  click  “Read  Data”  to  record  the  current  conditions. 

Then,  decrease  the  percent  of  gas  fuel,  or  make  a small  change  to  either  the  strain  rate  or  the  equivalence 
ratio.  Once  again,  small  changes  of  around  1 s'l  for  the  strain  rate,  and  0.005  for  the  equivalence  ratio 
are  necessary  to  prevent  instabilities.  When  the  new  conditions  have  been  selected,  click  on  the  “Go” 
arrow  to  adjust  to  these  settings. 

18.  During  testing,  it  is  important  to  keep  watch  of  the  current  burner  temperature  read-outs  on  the  screen 
and  to  make  adjustment  in  order  to  maintain  the  desired  burner  temperature.  Keep  the  upper  and  lower 
burner  temperatures  at  100°C  ± 5°C  by  adjusting  the  flow  rate  of  the  cooling/heating  water  into  the 
burner  and,  if  necessary,  the  variac  power  settings. 

19.  Repeat  steps  17  (or  16)  and  18,  changing  the  gas  composition  to  approach  the  extinction  limit.  If  the 
liquid  is  flammable  by  itself,  the  percent  of  gas  fuel  will  reduce  to  0 % and  the  LFL  tests  may  then  be 
conducted  by  reducing  the  equivalence  ratio  until  extinction  occurs.  For  R-245ca,  the  gas  fuel  cannot  be 
completely  eliminated,  and  the  maximum  R-245ca  vapor  amount  is  found  for  a particular  stretch  rate  and 
equivalence  ratio. 

20.  When  flame  has  been  extinguished,  press  the  “Read  Flows”  button  one  more  time  to  record  the  data 
at  the  extinction  condition.  Then  set  the  flow  control  button  to  “Off’,  and  click  on  the  “Go”  arrow  to 
stop  the  gas  flows.  If  testing  with  R-245ca,  shut  off  the  pump  and  close  the  valves  after  the  syringes. 

21.  Let  the  scrubber  continue  to  operate  for  at  least  five  minutes  after  the  flows  have  stopped,  in  order 
to  clean  any  combustion  gases  that  may  still  be  present  in  the  piping  system.  Check  the  scrubber 
drainage  container.  If  it  is  full,  or  close  to  being  full,  replace  with  an  empty  container  and  neutralize  the 
collected  liquid. 

22.  If  more  R-245ca  tests  are  to  be  run,  check  the  syringes  to  see  that  they  contain  enough  liquid.  Then 
repeat  steps  1 1 through  21.  Once  all  tests  are  completed,  the  system  may  be  shut  down. 


Shut  Down  Procedures: 

1.  Turn  off  valves  at  gas  bottles. 

2.  Empty  any  remaining  liquid  fuel  from  the  syringes. 
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3.  Depressurize  the  gas  lines  by  setting  all  MFCs  to  typical  flows  and  allowing  any  gas  left  in  the  lines 
to  flow  out. 

4.  When  lines  are  empty,  shut  off  valves  after  the  MFCs. 

5.  Set  the  MFCs  to  “off”  at  the  digital  control  box. 

6.  Turn  off  the  variacs  for  the  boiler,  heat  exchanger,  heated  lines  and  the  burner  heating  tape.  DO  NOT 
turn  off  heat  exchanger  water  flow  until  all  tubing  lines  have  cooled  completely. 
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APPENDIX  B.  Estimates  of  Test  Time  and  Material 


To  complete  a set  of  twelve  lean  flammability  tests  at  100  °C  and  with  humid  air  takes 
approximately  six  hours.  This  includes  two  hours  of  flow  calibration  time,  and  an  additional  four  hours 
of  testing  time.  For  a substance  similar  to  R-32,  about  450  grams  of  refrigerant  are  needed  to  complete 
the  testing.  These  estimates  assume  that  an  approximate  ignition  condition  is  known. 

Finding  the  critical  flammability  ratio  takes  about  seven  hours  when  using  a base  fuel  with 
known  characteristics.  This  allows  for  two  hours  to  calibrate  the  refrigerant  flow,  and  five  hours  to 
perform  the  set  of  twelve  tests.  For  an  inhibitor  like  R-125,  at  least  200  grams  are  necessary  to  complete 
the  tests.  Also,  about  400  grams  of  the  base  fuel  are  needed.  If  the  base  fuel  has  not  been  previously 
calibrated,  another  50  grams  of  fuel  and  two  hours  are  needed  to  perform  the  fuel  flow  calibration. 

These  estimates  are  for  one  set  of  tests  at  a predetermined  equivalence  ratio. 

The  liquid  tests  are  more  involved  and  require  additional  preparation  and  testing  time. 

Individual  tests  may  take  longer  than  30  minutes  each,  and  a set  of  experiments  takes  about  three  days  to 
complete,  depending  on  the  number  of  tests  needed  to  map  out  the  curve.  At  least  600  grams  of  liquid 
fuel,  and  1000  grams  of  vapor  base  fuel  are  needed  for  one  set  of  experiments,  assuming  the  calibrations 
have  already  been  completed.  This  estimate  is  based  on  testing  a slightly  flammable  liquid,  such  as 
R-245ca,  which  needs  to  be  mixed  with  a more  flammable  vapor  fuel  in  order  to  ignite  in  the  opposed- 
flow  burner. 


69 


APPENDIX  C:  Data  Acquisition  Program 

а.  Overview  of  Program 

The  refrigerant  experiments  were  conducted  using  TWINFLAM.VI,  a virtual  instrument  created  in 
National  Instruments  Lab  VIEW  3.1.1.  This  program  was  written  specifically  for  these  tests  and  uses  an 
IEEE  Standard  488  GPIB  interface  to  communicate  with  the  mass  flow  controllers.  A sample  print-out 
of  the  front  panel  screen  is  included  later  in  the  appendix.  Although  these  tests  were  conducted  using 
Lab  VIEW,  any  programming  tool  capable  of  reading  inputs  and  controlling  the  mass  flow  controllers 
(MFCs)  could  be  used.  A complete  copy  of  the  custom  program  is  available  from  the  authors.  The 
structure  of  the  program  is  outlined  in  the  following  steps: 

1 . The  computer  reads  the  inputs  for  the  types  of  gases  used,  and  for  the  desired  stretch  rate, 
equivalence  ratio,  and  percent  inhibitor,  and  calculates  the  respective  output  flow  rates  for  each  gas 
using  the  basic  equations. 

2.  The  computer  reads  the  inputs  for  the  relative  humidity,  barometric  pressure  and  the  current 
temperatures  in  the  lab  and  at  each  burner  outlet.  It  then  computes  the  settings  necessary  to  obtain 
the  correct  gas  flow  rates  for  the  conditions  chosen.  It  also  calculates  the  settings  for  the  nitrogen 
rotameters  based  on  their  calibration,  and  displays  these  so  that  they  may  be  set  manually  if  a 
nitrogen  co-flow  is  desired. 

3.  A digital  signal  is  sent  to  the  MFCs  to  begin  the  flows. 

4.  The  computer  then  monitors  the  actual  flow  rates  of  the  gases.  The  computer  queries  the  settings  at 
the  MFCs.  It  then  obtains  the  current  temperatures,  humidity,  and  barometric  pressure,  and 
computes  the  volumetric  flow  rate  of  each  gas  as  described  in  the  section  “Determination  of  Flow”. 
This  information  is  reported  on  the  screen. 

5.  Using  the  actual  flow  rates  of  the  gases,  current  values  for  the  strain  rate,  equivalence  ratio,  percent 
inhibitor,  and  the  velocity  at  the  burner  outlet  are  calculated.  This  process  of  reading  flow  rates  and 
calculating  current  conditions  is  repeated  every  few  seconds  and  the  screen  is  continually  updated 
accordingly. 

б.  When  the  “Read  Flows”  command  is  activated,  the  computer  stops  updating  conditions,  and  records 
all  the  current  data  to  the  chosen  file.  The  computer  then  goes  into  stand-by  mode  as  it  waits  for  the 
next  command.  This  means  that  the  gases  continue  to  flow  at  their  present  settings,  but  the  screen  is 
no  longer  updated.  The  computer  waits  for  the  user  to  either  indicate  new  inputs  and  run  the 
program  again,  or  to  shut  off  gas  flows  and  end  the  program. 
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b.  Determination  of  Flow 


The  mass  flow  controllers  set  and  measure  flow  as  an  absolute  amount  of  mass  that  passes  through  the 
controller  in  a given  amount  of  time.  This  mass  flow  rate  measured  at  the  controllers  is  changed  to  the 
volumetric  flow  rate  through  the  burner  by  the  following  series  of  conversions: 

1 . The  mass  flow  controller  measures  the  flowing  gas  in  units  of  mass  per  unit  of  time.  Internally,  the 
MFC  converts  this  measurement  into  a volume  flow  rate  for  nitrogen  (the  default  gas)  at  standard 
temperature  and  pressure.  This  number  is  then  transmitted  to  the  computer  via  a GPIB  connection. 

2.  Next,  the  computer  program  changes  the  nitrogen  flow  rate  into  the  corresponding  flow  rate  for  the 
actual  gas  being  used.  Prior  to  conducting  the  experiments,  each  flow  meter  is  calibrated  with  its 
designated  gas  at  standard  temperature  and  pressure,  using  a bubble  flow  meter.  This  calibration 
determines  the  linear  relationship  between  the  default  nitrogen  flow  rate,  and  the  measured  flow  of 
the  actual  gas,  as  well  as  any  zero  offset  in  the  meter  itself.  The  computer  program  uses  these  values 
to  calculate  the  flow  rate  of  the  actual  gas  at  standard  temperature  and  pressure  using  the  equation: 

Qstp  = m x On  2 + b 

where  m is  about  equal  to  the  ratio  of  the  nitrogen  specific  heat  times  density  to  the  actual  gas 
specific  heat  times  density,  and  b is  the  zero  offset. 

3.  The  flow  rate  for  the  gas  must  then  be  corrected  from  standard  temperature  and  pressure,  to  the 
measured  pressure  in  the  laboratory  and  the  temperature  at  the  burner  outlet.  This  is  done  using  the 
ideal  gas  law. 

( Tactual']  f PSTP 

= MTrJte 

4.  Finally,  the  flow  is  corrected  for  humidity  effects  using  the  saturation  vapor  pressure  of  the  water. 
The  saturation  vapor  pressure  is  calculated  as: 

Py  = (-16.558)  +(2.85357)  - (0.079649J2)  + (0.001286T3) 

where  T is  temperature  measured  in  degrees  Kelvin  and  Pv  is  in  torr.  The  final  corrected  flow  rate  is 
computed  by: 


Qcorr  — Q 


P actual 


P actual  — (Pv)(RH)J 


where  RH  is  the  relative  humidity  and  P actual  is  the  measured  barometric  pressure  in  torr. 
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c.  Basic  Equations 

The  following  equations  describe  the  relationships  between  the  quantities  measured  in  these  experiments. 
The  program  uses  these  equations  to  calculate  the  initial  flow  rates  to  be  set  (Qrj,  Qr2,  and  Qa[r)  when 
the  desired  equivalence  ratio,  stretch  rate  and  percent  inhibitor  are  given  as  inputs.  It  then  applies  these 
equations  to  calculate  the  resulting  equivalence  ratio,  stretch  rate  and  percent  inhibitor  based  on  the 
measured,  updated  flow  rates. 


R1  = more  flammable  refrigerant 
R2  = less  flammable  refrigerant 
K = stoichiometric  air/fuel  ratio 


Velocity: 


( 


Stretch  Rate:  SR  = 


v 


Percent  Inhibitor:  % Inhibitor  = 
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d.  Front  Panel 


The  front  panel  of  the  program  TWINFLAM.VI  is  shown  on  the  following  page.  It  is  composed  of  the 
elements  listed  below.  The  front  panel  is  the  only  program  screen  used  when  conducting  the 
flammability  tests  and  contains  all  the  computer  controls  needed.  All  inputs,  and  the  resulting  calculated 
outputs  appear  on  this  screen,  as  well. 


Manual  Inputs: 

File  Name 

Fuel  Selector  (choose  gases  used  during  testing) 

Run  Number  (differentiate  between  tests  in  same  file) 

Stretch  Rate  (s~  1 ) 

Equivalence  Ratio 

Percent  Inhibitor  (volume  percent  of  fuel) 

Relative  Humidity  (percent) 

Nitrogen  Co-flow 

Automatic  Inputs: 

Barometric  Pressure  (torr) 

Upper  Burner  Temperature  (°C) 

Lower  Burner  Temperature  (°C) 

Room  Temperature  (°C) 

Controls: 

Go  (runs  program) 

Read  Flows  (records  data  to  file) 

All  MFCs  On/Off 

Individual  MFCs  On/Off 

Calculated  Outputs: 

Actual  Equivalence  Ratio 

Actual  Stretch  Rate 

Percent  Inhibitor 

Total  Fuel  Flow  (LPM) 

Total  Air  Flow  (LPM) 

Total  Inhibitor  Flow  (LPM) 

Velocity  of  Flow  (cm/s) 

Temperature  Difference  (between  upper  and  lower  burner  outlets) 

Nitrogen  Rotameter  settings 

Nitrogen  Flow  (LPM) 
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TWINFLAM.VI 
03/27/98  05:47  PM 


GO 

ON 

o 


MFC  Controls 

o c c c o o 


File  name 
°0  c:\dalal  .doc 


31 


Run  U 


FUEL  SELECTOR 

S : R32 /R12  5 /Air  T7  2 


Strain  Rate 


5170.0 


Equiv.  Ratio 


Nitrogen  Fraction 

^jOTJcF 


% Inhibitor 


0.0 


RH 


% |49.80 


N2  Meter  Flow  (L/m) 


0.000 


Top  Burner 
Rotameter 

1;  66.00  ll 

Bottom  Burner 
Rotameter 


62.50 


Total  Fuel  (L/m) 


0.0000 


Total  Air  (L/m) 


0.000 


Total  Inhibitor  (L/m) 


0.0000 


Velocity  (cni/s) 

il  0.00  II 


Temperature  (C) 


23.89 


Temperatures 


Equiv.  Ratio 


1.0002  I <D 


Strain  Rate 


69.98 


Inhibitor  % 


0.00 


Pressure  (Torr) 


749.997  BP 


Delta  T 


Read 

Flows 


Top  Burner  (Q 


100.1 


Lower  Burner  (Q 


99.8 
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APPENDIX  D.  C-H-O-F  Chemical  Kinetics  Mechanism  [7,  16,  17,  18] 


CHEMKIN  INTERPRETER  OUTPUT:  CHEMKIN-II  Version 
3.9  Aug.  1994 

temperature  element  count 


species 

phase 

MW 

low 

high 

HOCNFAr 

1.02 

G 0 

31.99880 

200 

3500 

0 

2 0 

0 

0 

0 

2.  C02 

G 0 

44.00995 

200 

3500 

0 

2 1 

0 

0 

0 

3.  H20 

G 0 

18.01534 

200 

3500 

2 

1 0 

0 

0 

0 

4.  CO 

G 0 

28.01055 

200 

3500 

0 

1 1 

0 

0 

0 

5.CH4 

G 0 

16.04303 

200 

3500 

4 

0 1 

0 

0 

0 

6.  CH2F2 

G 0 

52.02389 

300 

3000 

2 

0 1 

0 

2 

0 

7.  H2 

G 0 

2.01594 

200 

3500 

2 

0 0 

0 

0 

0 

8.  OH 

G 0 

17.00737 

200 

3500 

1 

1 0 

0 

0 

0 

9.  H 

G 0 

1.00797 

200 

3500 

1 

0 0 

0 

0 

0 

10.0 

G 0 

15.99940 

200 

3500 

0 

1 0 

0 

0 

0 

11.  HF 

G 0 

20.00637 

300 

5000 

1 

0 0 

0 

1 

0 

12.  CF2:0 

G 0 

66.00735 

300 

3000 

0 

1 1 

0 

2 

0 

13.  CF4 

G 0 

88.00475 

300 

3000 

0 

0 1 

0 

4 

0 

14.CH3 

G 0 

15.03506 

200 

3500 

3 

0 1 

0 

0 

0 

15.CH2 

G 0 

14.02709 

200 

3500 

2 

0 1 

0 

0 

0 

16  . CH 

G 0 

13  01912 

200 

3500 

1 

0 1 

0 

0 

0 

17.  H02 

G 0 

33.00677 

200 

3500 

1 

2 0 

0 

0 

0 

18.  H202 

G 0 

34.01474 

200 

3500 

2 

2 0 

0 

0 

0 

19.  CH2(S) 

G 0 

14.02709 

200 

3500 

2 

0 1 

0 

0 

0 

20.  C 

G 0 

12.01115 

200 

3500 

0 

0 1 

0 

0 

0 

21.C2H6 

G 0 

30.07012 

200 

3500 

6 

0 2 

0 

0 

0 

22.  C2H5 

G 0 

29.06215 

200 

3500 

5 

0 2 

0 

0 

0 

23.  C2H4 

G 0 

28.05418 

200 

3500 

4 

0 2 

0 

0 

0 

24.  C2H3 

G 0 

27.04621 

200 

3500 

3 

0 2 

0 

0 

0 

25.  C2H2 

G 0 

26.03824 

200 

3500 

2 

0 2 

0 

0 

0 

26.  C2H 

G 0 

25.03027 

200 

3500 

1 

0 2 

0 

0 

0 

27.  CH30H 

G 0 

32.04243 

200 

3500 

4 

1 1 

0 

0 

0 

28.  CH30 

G 0 

31.03446 

300 

3000 

3 

1 1 

0 

0 

0 

29.  CH20H 

G 0 

31.03446 

200 

3500 

3 

1 1 

0 

0 

0 

30.  CH20 

G 0 

30.02649 

200 

3500 

2 

1 1 

0 

0 

0 

31.  HCO 

G 0 

29.01852 
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1 

1 1 

0 

0 

0 

32.  CH2CO 

G 0 

42.03764 

200 

3500 

2 

1 2 

0 

0 

0 

33.  HCCO 

G 0 

41.02967 

300 

4000 

1 

1 2 

0 

0 

0 

34.  HCCOH 

G 0 

42.03764 

300 

5000 

2 

1 2 

0 

0 

0 

35.  F 

G 0 

18.99840 

300 

5000 

0 

0 0 

0 

1 

0 

36.  CH3F 

G 0 

34.03346 

300 

3000 

3 

0 1 

0 

1 

0 

37.  CHF3 

G 0 

70.01432 

300 

3000 

1 

0 1 

0 

3 

0 

38.  CH2F 

G 0 

33.02549 

300 

3000 

2 

0 1 

0 

1 

0 

39.  CHF2 

G 0 

51.01592 

300 

3000 

1 

0 1 

0 

2 

0 

40.  CF3 

G 0 

69.00635 

300 

3000 

0 

0 1 

0 

3 

0 

41.CHF 

G 0 

32.01752 

300 

3000 

1 

0 1 

0 

1 

0 

42.  CF2 

G 0 

50.00795 

300 

3000 

0 

0 1 

0 

2 

0 

43.  CF 

G 0 

31.00955 

300 

3000 

0 

0 1 

0 

1 

0 

44.  CF30 

G 0 

85.00575 

300 

3000 

0 

1 1 

0 

3 

0 

45.  CHF:0 

G 0 

48.01692 

300 

3000 

1 

1 1 

0 

1 

0 

46.  CF:0 

G 0 

47.00895 

300 

3000 

0 

1 1 

0 

1 

0 

47.  CHFCO 

G 0 

60.02807 

300 

3000 

1 

1 2 

0 

1 

0 

48.  CF2CO 

G 0 

78.01850 

300 

3000 

0 

1 2 

0 

2 

0 

49.  FCCO-E 

G 0 

59.02010 

300 

3000 

0 

1 2 

0 

1 

0 

50.  CH3-CH2F 

G 0 

48.06055 

300 

1500 

5 

0 2 

0 

1 

0 

51.  CH3-CHF2 

G 0 

66.05098 

300 

1500 

4 

0 2 

0 

2 

0 

52.  CH3-CF3 

G 0 

84.04141 

300 

1500 

3 

0 2 

0 

3 

0 

53.  CH2F-CH2F 

G 0 

66.05098 

300 

1500 

4 

0 2 

0 

2 

0 

54.  CH2F-CHF2 

G 0 

84.04141 

300 

3000 

3 

0 2 

0 

3 

0 

55.  CH2F-CF3 

G 0 

102.03184 

300 

1500 

2 

0 2 

0 4 

0 

56.  CHF2-CHF2 

G 0 102.03184 

300 

3000 

2 

0 2 

0 4 

0 

57.  CHF2-CF3 

G 0 

120.02227 

300 
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1 

0 2 

0 

5 

0 

58.  CF3-CF3 

G 0 

138.01270 

300 

3000 

0 

0 

2 

0 

6 

0 

59.  CH3-CHF 

G 0 

47.05258 

300 
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4 

0 

2 

0 

1 

0 

60.  CH3-CF2 

G 0 

65.04301 
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3 

0 

2 

0 

2 
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G 0 
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0 
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0 
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G 0 
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0 

2 

0 

2 
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G 0 
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0 

2 

0 
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G 0 
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2 
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2 
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2 

0 

2 
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3 
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1 

0 

2 
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0 
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0 
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0 
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0 

2 

0 
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1 

0 

2 

0 

3 
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G 0 

100.01590 

300 

3000 

0 
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0 
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G 0 
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2 

0 

2 

0 

1 
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G 0 
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2 

0 

2 

0 

1 
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G 0 
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1 

0 

2 

0 

2 
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G 0 
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1 

0 

2 

0 

2 
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G 0 
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0 

0 

2 

0 
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G 0 
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1 

0 

2 

0 
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G 0 
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0 

0 

2 

0 

2 
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G 0 

28.01340 

300 

5000 

0 

0 

0 

2 

0 

0 

83.  AR 

G 0 

39.94800 

300 
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0 

0 

0 

0 

0 

1 
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APPENDIX  E.  Detailed  Drawings  of  Industrial  Burner  (dimensions  in  inches  and  mm) 
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APPENDIX  F:  Drawings  of  Research  Burner 
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